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(54) MR imaging of synchronous spin motion and strain waves 



(57) A scan using an NMR imaging system is carried 
out while applying an oscillating stress (319) to the object 
being imaged. An alternating magnetic field gradient 
(315) synchronized with the applied stress (319) is 
employed in the NMR imaging pulse sequence (A) to 

FIG. 3 

315 v 



detect and measure synchronous spin motion (309A) 
throughout the field of view. The direction of the alternat- 
ing gradient (315) and/or the applied stress (319) may 
be changed to measure and image the elastic properties 
of the object. 
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Description 

Background of the Invention 

5 The field of the invention is nuclear magnetic resonance imaging methods and systems. More particularly, the inven- 

tion relates to the enhancement of MR image contrast. 

The physician has many diagnostic tools at his or her disposal which enable detection and localization of diseased 
tissues. These include x-ray systems that measure and produce images indicative of the x-ray attenuation of the tissues 
and ultrasound systems that detect and produce images indicative of tissue echogenicity and the boundaries between 

10 structures of differing acoustic properties. Nuclear medicine produces images indicative of those tissues which absorb 
tracers injected into the patient, as do PET scanners and SPECT scanners. And finally, magnetic resonance imaging 
("MRI") systems produce images indicative of the magnetic properties of tissues. It is fortuitous that many diseased 
tissues are detected by the physical properties measured by these imaging modalities, but it should not be surprising 
that many diseases go undetected. 

is Historically, one of the physician's most valuable diagnostic tools is palpation. By palpating the patient a physician 
can feel differences in the compliance of tissues and detect the presence of tumors and other tissue abnormalities. 
Unfortunately, this valuable diagnostic tool is limited to those tissues and organs which the physician can feel, and many 
diseased internal organs go undiagnosed unless the disease happens to be detectable by one of the above imaging 
modalities. Tumors (e.g. of the liver) that are undetected by existing imaging modalities and cannot be reached for 

20 palpation through the patient's skin and musculature, are often detected by surgeons by direct palpation of the exposed 
organs at the time of surgery. Palpation is the most common means of detecting tumors of the prostate gland and the 
breast, but unfortunately, deeper portions of these structures are not accessible for such evaluation. An imaging system 
that extends the physician's ability to detect differences in tissue compliance throughout a patient's body would extend 
this valuable diagnostic tool. 

25 Any nucleus which possesses a magnetic moment attempts to align itself with the direction of the magnetic field in 
which it is located. In doing so, however, the nucleus precesses around this direction at a characteristic angular frequency 
(Larmor frequency) which is dependent on the strength of the magnetic field and on the properties of the specific nuclear 
species (the magnetogyric constant y of the nucleus). Nuclei which exhibit this phenomena are referred to herein as 
"spins", and materials which contain such nuclei are referred to herein as "gyromagnetic". 

30 When a substance such as human tissue is subjected to a uniform magnetic field (polarizing field B 0 ), the individual 
magnetic moments of the spins in the tissue attempt to align with this polarizing field, but precess about it in random 
order at their characteristic Larmor frequency. A net magnetic moment M z is produced in the direction of the polarizing 
field, but the randomly oriented magnetic components in the perpendicular, or transverse, plane (x-y plane) cancel one 
another. If, however, the substance, or tissue, is subjected to a magnetic field (excitation field Bi) which is in the x-y 

35 plane and which is near the Larmor frequency, the net aligned moment, Mz, may be rotated, or "tipped", into the x-y 
plane to produce a net transverse magnetic moment M t , which is rotating, or spinning, in the xy plane at the Larmor 
frequency. The practical value of this phenomenon resides in the signal which is emitted by the excited spins after the 
excitation signal is terminated. There are a wide variety of measurement sequences in which this nuclear magnetic 
resonance ("NMR") phenomena is exploited. 

40 When utilizing NMR to produce images, a technique is employed to obtain NMR signals from specific locations in 
the subject. Typically, the region which is to be imaged (region of interest) is scanned by a sequence of NMR measure- 
ment cycles which vary according to the particular localization method being used. The resulting set of received NMR 
signals are digitized and processed to reconstruct the image using one of many well known reconstruction techniques. 
To perform such a scan, it is, of course, necessary to elicit NMR signals from specific locations in the subject. This is 

45 accomplished by employing magnetic fields (G x , G y and G z ) which are superimposed on the polarizing field B 0 , but 
which have a gradient along the respective x, y and z axes. By controlling the strength of these gradients during each 
NMR cycle, the spatial distribution of spin excitation can be controlled and the location of the resulting NMR signals can 
be identified. 

It is well known that NMR can be used to detect and image the movement of spins. As disclosed in U.S. patent No. 
so Re. 32,701 entitled "NMR Scanner With Motion Zeugmatographf acquired NMR signals can be sensitized to detect 
moving spins by applying a bipolar magnetic field gradient at the proper moment in each NMR measurement sequence. 
The phase of the resulting NMR signal measures the velocity of spins along the direction of the motion sensitizing 
magnetic field gradient. With more complex motion sensitizing magnetic field gradients, higher orders of motion, such 
as acceleration and jerk can also be measured with this method. 

55 

Summary of the Invention 

The present invention relates to the imaging of mechanical characteristics of an object, and more particularly, the 
production of an NMR image in which the intensity pattern reveals the mechanical properties of the object and in some 
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applications, the pattern of propagation of mechanical strain waves within it. The invention includes the acquisition of 
NMR imaging data using an NMR pulse sequence in which a magnetic field gradient is alternated in polarity to phase 
encode moving spins, and a stimulus is applied to the object being imaged and is oscillated in magnitude in synchronism 
with the alternating magnetic field gradient. A phase image is reconstructed from the acquired NMR imaging data which 

5 depicts the synchronous movement of spins throughout the object resulting from the applied stimulus. Mechanical char- 
acteristics of the object are revealed by analysis of such images. 

A general object of the invention is to produce an image which reveals the mechanical properties of an object. The 
applied stimulus may be a mechanical disturbance which sets up a traveling wave in the gyromagnetic medium that 
causes cyclic strain, or displacements of the spins. By synchronizing the alternating magnetic field gradient with the 

10 oscillatory applied stress, the movement of the spins in relation to the applied force is revealed by the phase of the 
acquired NMR signals. At locations where the spins move in phase with the applied stress the phase signal is maximum, 
at one polarity, and at locations where the spins oscillate out of phase with the applied stress the phase signals are 
maximum at the opposite polarity. These phase signals may be employed to modulate pixel intensity in a reconstructed 
image and the strain due to the applied stress throughout the imaged object is seen as a wave pattern. Changes in this 

15 wave pattern reveal the boundaries of regions of differing mechanical properties within the object. 

Another object of the invention is to produce a series of images which reveal the strain wave propagation through 
an object. By synchronizing the application of the oscillating applied stress with the alternating magnetic field gradient 
at a series of different phase relationships, a corresponding series of phase images may be reconstructed showing the 
movement of the strain wave pattern through the object. 

20 Yet another object of the invention is to produce an image which indicates the strain in an object produced by an 
applied force. A first phase image is produced with the oscillating stress applied, and a second phase image is produced 
without the oscillating stress applied. The difference between these two phase images is calculated and provides an 
image indicative of the strain in the object. These strain images provide a measure of a fundamental quantity of the 
medium in the form of displacement amplitude as a strain wave propagates through it. As a result, it becomes possible 

25 to calculate a number of fundamental mechanical properties pertaining to the viscoelastic nature of the gyromagnetic 
medium. They include, for example, velocity of propagation, Young's modulus and other moduli. 

Another object of the invention is to produce images of the object indicative of its viscoelastic properties. By producing 
two phase images with the oscillating stress synchronized with the alternating magnetic field gradient at two different 
phase relationships, the wavelength of the propagated strain can be calculated at each image pixel. From this, the 

30 velocity of propagation can be calculated at each image pixel as well as Young's modulus. 

Another object of the invention is to produce an image which indicates Poisson's ratio throughout the object. A first 
phase image is produced with the alternating gradient applied in the same direction as the oscillating stress, and a 
subsequent phase image is produced with the alternating gradient applied in an orthogonal direction. Poisson's ratio 
may be estimated at each pixel location from the resulting phase images. Since the velocity and wavelength of the 

35 propagated strain may also be measured, the viscous properties of the gyromagnetic medium can also be estimated. 
Using Young's modulus and Poisson's ratio, other moduli such as the shear modulus and the bulk modulus can be 
calculated and imaged. 

Another general object of the invention is to produce an image of a patient which indicates the compliance of tissues. 
By acquiring a plurality of phase images of the patient with different arrangements of the oscillatory stress and the 
40 alternating magnetic field gradient, an image indicative of tissue compliance can be produced. This image provides the 
clinician with information comparable to that obtained by manually palpating the tissues. 

Yet another object of the invention is to produce an NMR image in which contrast is enhanced for spins having 
oscillatory motion. By synchronizing the alternating magnetic field gradient with oscillatory movement of spins in the 
object being imaged, a phase image may be reconstructed in which the signals from the oscillatory spins are enhanced 
45 relative to stationary spins and randomly moving spins. 

The foregoing and other objects and advantages of the invention will appear from the following description. In the 
description, reference is made to the accompanying drawings which form a part hereof, and in which there is shown by 
way of illustration a preferred embodiment of the invention. Such embodiment does not necessarily represent the full 
scope of the invention, however, and reference is made therefore to the claims herein for interpreting the scope of the 
so invention. 

Brief Description Of The Drawings 

Fig. 1 is a block diagram of an NMR system which employs the present invention; 
55 Fig. 2 is an electrical block diagram of the transceiver which forms part of the NMR system of Fig. 1 ; 

Fig. 3 is a graphic representation of a pulse sequence performed by the NMR system of Fig. 1 to practice the 
preferred embodiment of the invention; 
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Fig. 4 is a flow chart which illustrates how the NMR system of Fig. 1 reconstructs a series of strain images in 
accordance with a first preferred embodiment of the invention using NMR data acquired with the pulse sequence 
of Fig. 3; 

Fig. 5 is a flow chart which illustrates how the NMR system of Fig. 1 produces images indicative of the elastic 
properties in accordance with a second preferred embodiment of the invention using NMR data acquired with the 
pulse sequence of Fig. 3; 

Figs. 6a-c are flow charts of three different methods used to produce phase difference images in the methods of 
Figs. 4 and 5; 

Figs. 7a and 7b are elevation and top views respectively of a transducer that is used in the NMR system of Fig. 1 
to practice the present invention; 

Fig. 8 is a graphic illustration of the phase accumulation in an NMR signal acquired with the pulse sequence of Fig. 3; 
Fig. 9 is a graphic illustration of phase accumulation with the alternating gradient waveform reversed 180°; 
Figs. 10a and 10b are graphic illustrations of alternative methods for applying the oscillating stress relative to the 
alternating gradient waveform in the pulse sequence of Fig. 3; and 

Fig. 1 1 is a pictorial representation of a human leg illustrating an alternative method for stimulating synchronous 
spin motion in blood vessels. 

General Description Of The Invention 

In general, the physical properties of tissue are measured by applying a stress (e.g. tension, pressure, or shear) 
and observing the resulting strain (e.g. elongation, compression, rotation). By measuring the resulting strain, elastic 
properties of the tissue such as Young's modulus, Poisson's ratio, the shear modulus, and the bulk modulus, can be 
calculated. By applying the stress in all three dimensions and measuring the resulting strain, the elastic properties of 
the tissue can be completely defined. 

By observing the rate at which the strain decreases as a function of distance from the stress producing source, the 
attenuation of the strain wave can be estimated. From this, the viscous properties of the gyromagnetic medium may be 
estimated. The dispersion characteristics of the medium can be estimated by observing the speed and attenuation of 
the strain waves as a function of their frequency. Dispersion is potentially a very important parameter for characterizing 
tissues in medical imaging applications. 

The present invention provides a means for measuring the strain in gyromagnetic materials such as tissues using 
NMR methods and apparatus. The Larmor equation is given by 

co = yg (1) 

where co is the angular resonant frequency of the NMR signal produced by spins having a characteristic gyromagnetic 
ratio of y when placed in a magnetic field having a density and direction B. The magnetic field vector ~B can be broken 
down into two components 

B=B 0 + rG r (2) 

where B 0 is the polarizing magnetic field, r is the location of the spins, and G r is the magnetic field gradient. 

Since the angular frequency of the NMR signals produced by the spins is the rate of change of their phase, the 
phase of the spin signals as a function of time is as follows: 

♦(Q-Ja>(0c#. (3) 

Substituting equation (1), the relationship between NMR signal phase and the applied gradient field is obtained 

Mt) = YiG(t)r(f)dt\ (4) 

where G and r are expressed as functions of time (f) for obtaining a general expression of the Larmor equation. This 
equation indicates that the NMR signal produced by moving spins will accrue a phase shift relative to that accrued by 
static spins when in the presence of a magnetic field gradient. 

If an oscillating stress is applied to tissue along the direction r at an angular frequency cop, a wave is launched and 
spins are displaced by amounts determined by the elastic properties of the tissue. If it is assumed that this propagation 
occurs without loss, the displacement (A) of spins at location (r) may be expressed as follows: 



A = A 0 cos (co p t + kr +e), 



(5) 
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where A 0 is the maximum displacement produced by the applied stress, k is the wave number, and e is the phase offset 
of the spin displacement relative to the applied oscillating stress. The wave number k is equal to 2n radians divided by 
the wavelength (X) of the propagated wave, and if it is assumed that the spin displacement occurs for just one cycle (t=0 
to T) of the applied oscillating stress, then the NMR signal produced by the spin will accumulate a phase indicated by 
5 the following expression: 



t=T 



10 



40 



<K0=yJ G[f)A Q cos[<o p t+kr+e)df (6) 



f=0 



If the magnetic field gradient G(t') is constant during this time period, no phase signal will be accumulated. However, if 
the magnetic field gradient G(t') is synchronized with the applied stress and is switched in polarity half way through the 
time period (T), the phase of the NMR signal (<|)) at the completion of the time period will be proportional to the displace- 
ment of the spins. This displacement along the r direction is the strain which results from the applied stress along the 
15 same direction r. 

The phase accumulation of the NMR signal due to spins in synchronous motion is illustrated in Fig. 8. An alternating 
magnetic field gradient that completes one bi-polar cycle during the time period T is indicated by waveform 1 0. Spins in 
synchronous motion and perfectly in phase with the alternating magnetic field gradient 1 0 is illustrated by the waveform 
12, and the resulting accumulated phase in the NMR signal produced by these spins is indicated by the graph 14. On 

20 the other hand, synchronous spin motion 180° out of phase with the alternating gradient 10 is indicated by waveform 
13 and the resulting phase accumulation is shown in graph 15. In contrast, Fig. 9 illustrates at graph 16 the phase of 
NMR signals produced by stationary spins, and graph 18 indicates the accumulated phase of constantly moving spins. 
At the completion of one cycle of the alternating gradient 10, the accumulated phase of stationary spins is nulled and 
the accumulated phase of constantly moving spins is proportional to their velocity along the direction of the gradient. 

25 Because the displacements are very small, the gradient field G(t') is usually oscillated in synchronism with the 
applied stress for several cycles before the NMR signal is acquired. This enables the accumulated phase (()>) to reach 
a significant amount. A second NMR signal is then acquired with an identical pulse sequence, but without the application 
of the oscillatory stress. This second signal provides a reference phase (tyfd which may be subtracted from the accu- 
mulated phase <|> A to yield a value indicative of longitudinal strain (S|_). 

30 These NMR measurements can be made, of course, with imaging gradients applied, and a strain image indicative 
of the longitudinal strain at each pixel may be reconstructed. This strain image has pixel values as follows: 

S L (t,r) = A 0 cos(co p t + kr). (7) 

35 If the measurement is repeated with the oscillating stress in the same direction, but with the alternating gradient 
G(t') oriented along each of the orthogonal axes, the displacements, along these axes can be determined and the orthog- 
onal strain (S T ) calculated. From this information Poisson's ratio (a) can be calculated as follows: 



a = S T /S L (8) 

Further information can be learned about the elastic properties of tissues by changing the phase relationship of 
the applied oscillatory stress and the synchronized alternating motion gradient. For example, if the phases are offset 
n/2 radians, a second strain image So is produced having the following pixel values: 

45 S 0 (t,r) = A 0 cos(co t + kr + n/2) 

(9) 

= Aosin(co p t + kr) 

The gradient in each of the strain images S L and So may be defined as follows: 

so vS L = - r t/+- r t; (10) 

x y 

55 x y 

Using these gradients, the wave number (k) of the propagated wave at each image pixel can be calculated as follows: 



k = VIvSJ ■ |vS L |+|vS 0 | ■ \vS 0 \/Js 2 0 + S 2 L (1 2) 
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Knowing the frequency (f) of the applied oscillatory stress, the propagation velocity (c) can then be calculated at each 
image pixel as long as the viscosity effects are not significant. 



10 



15 



20 



25 



30 



35 



X = 2n/k 

C = fX 



(13) 



If the density (p) of the gyromagnetic medium is known, the propagation velocity (c) can be used to calculate Young's 
modulus (Y): 



Y = c2p. 



(14) 



Since compliance is the inverse of Young's modulus (i.e.1/Y), the compliance of the gyromagnetic medium may also 
be calculated. An image in which pixel intensity is determined by the calculated compliance has diagnostic value in 
medicine because such an image displays what a physician feels when manually palpating tissue. With the knowledge 
of Poisson's ratio (a) and the Young's modulus (Y) all other values of moduli, namely the shear modulus Qui) and the bulk 
modulus (p) may be calculated, since only two of these four elastic properties are actually independent. Bulk modulus 
(P) can be expressed from Young's modulus and Poissons ratio, in the following way: 



Y = 3P(1 -2a). 



Shear Modulus (jli) can be expressed as follows: 



Y = 2jx(1 + a). 



(15) 



(16) 



It should be apparent that the oscillatory stress may be applied in three orthogonal directions and the synchronized 
gradient field may also be applied in three separate orthogonal directions for each applied stress direction. Spin dis- 
placements (A, r|, L) in all three directions may be measured, and all the components of the strain dyadic (second-rank 
tensor field), given in the following matrix may be measured: 



(17) 



(dA 
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dx 


dx 


dx 
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dy 
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dL 
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dz 


dz) 



With these measurements it is possible to characterize all the elastic properties of the gyromagnetic medium under 
40 investigation using the above-described calculations. 

Description Of The Preferred Embodiment 

Referring first to Fig. 1, there is shown the major components of a preferred NMR system which incorporates the 
45 present invention and which is sold by the General Electric Company under the trademark "SIGNA". The operation of 
the system is controlled from an operator console 100 which includes a console processor 101 that scans a keyboard 
102 and receives inputs from a human operator through a control panel 103 and a plasma display/touch screen 104. 
The console processor 101 communicates through a communications link 116 with an applications interface module 

1 1 7 in a separate computer system 1 07. Through the keyboard 1 02 and controls 1 03, an operator controls the production 
so and display of images by an image processor 1 06 in the computer system 1 07, which connects directly to a video display 

1 18 on the console 100 through a video cable 105. 

The computer system 107 includes a number of modules which communicate with each other through a backplane. 
In addition to the application interface 1 1 7 and the image processor 1 06, these include a CPU module 1 08 that controls 
the backplane, and an SCSI interface module 109 that connects the computer system 107 through a bus 1 10 to a set 
55 of peripheral devices, including disk storage 1 1 1 and tape drive 1 1 2. The computer system 1 07 also includes a memory 
module 113, known in the art as a frame buffer for storing image data arrays, and a serial interface module 114 that 
links the computer system 107 through a high speed serial link 115 to a system interface module 120 located in a 
separate system control cabinet 122. 
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The system control 122 includes a series of modules which are connected together by a common backplane 118. 
The backplane 1 18 is comprised of a number of bus structures, including a bus structure which is controlled by a CPU 
module 119. The serial interface module 120 connects this backplane 1 18 to the high speed serial link 1 15, and pulse 
generator module 121 connects the backplane 118 to the operator console 100 through a serial link 125. It is through 
5 this link 125 that the system control 122 receives commands from the operator which indicate the scan sequence that 
is to be performed. 

The pulse generator module 121 operates the system components to carry out the desired scan sequence. It pro- 
duces data which indicates the timing, strength and shape of the RF pulses which are to be produced, and the timing 
of and length of the data acquisition window. The pulse generator module 121 also connects through serial link 126 to 

10 a set of gradient amplifiers 1 27, and it conveys data thereto which indicates the timing and shape of the gradient pulses 
that are to be produced during the scan. 

In the preferred embodiment of the invention the pulse generator module 121 also produces sync pulses through 
a serial link 128 to a wave generator and amplifier 129. The wave generator produces a sinusoidal voltage which is 
synchronized to the frequency and phase of the received sync pulses and this waveform is output through a 50 watt, dc 

15 coupled audio amplifier. A frequency in the range of 20 Hz to 1000 Hz is produced depending on the particular object 
being imaged, and it is applied to a transducer 130. The transducer 130 will be described in more detail below, and its 
structure depends on the particular anatomy being measured and imaged. In general, however, the transducer 130 
produces a force, or pressure, which oscillates in phase with the sync pulses produced by the pulse generator module 
121 and creates an oscillating stress in the gyromagnetic media (i.e. tissues) to which it is applied. 

20 And finally, the pulse generator module 121 connects through a serial link 132 to scan room interface circuit 133 
which receives signals at inputs 135 from various sensors associated with the position and condition of the patient and 
the magnet system. It is also through the scan room interface circuit 133 that a patient positioning system 134 receives 
commands which move the patient cradle and transport the patient to the desired position for the scan. 

The gradient waveforms produced by the pulse generator module 121 are applied to a gradient amplifier system 

25 127 comprised of G x , G y and G z amplifiers 136, 137 and 138, respectively. Each amplifier 136, 137 and 138 is utilized 
to excite a corresponding gradient coil in an assembly generally designated 139. The gradient coil assembly 139 forms 
part of a magnet assembly 141 which includes a polarizing magnet 1 40 that produces either a 0.5 or a 1 .5 Tesla polarizing 
field that extends horizontally through a bore 142. The gradient coils 139 encircle the bore 142, and when energized, 
they generate magnetic fields in the same direction as the main polarizing magnetic field, but with gradients G x , G y and 

30 G z directed in the orthogonal x-, y- and z-axis directions of a Cartesian coordinate system. That is, if the magnetic field 
generated by the main magnet 140 is directed in the z direction and is termed B 0 , and the total magnetic field in the z 
direction is referred to as B z ,then G x =dB z /dx, G y =3B z /8y and G z =d B z /dz , and the magnetic field at any point (x,y,z) 
in the bore of the magnet assembly 141 is given by B(x,y,z)=B 0 +G x x+G y y+G z z . The gradient magnetic fields are 
utilized to encode spatial information into the NMR signals emanating from the patient being scanned, and as will be 

35 described in detail below, they are employed to measure the microscopic movement of spins caused by the pressure 
produced by the transducer 130. 

Located within the bore 142 is a circular cylindrical whole-body RF coil 152. This coil 152 produces a circularly 
polarized RF field in response to RF pulses provided by a transceiver module 150 in the system control cabinet 122. 
These pulses are amplified by an RF amplifier 1 51 and coupled to the RF coil 1 52 by a transmit/receive switch 1 54 which 

40 forms an integral part of the RF coil assembly. Waveforms and control signals are provided by the pulse generator module 
121 and utilized by the transceiver module 150 for RF carrier modulation and mode control. The resulting NMR signals 
radiated by the excited nuclei in the patient may be sensed by the same RF coil 152 and coupled through the trans- 
mit/receive switch 154 to a preamplifier 153. The amplified NMR signals are demodulated, filtered, and digitized in the 
receiver section of the transceiver 1 50. The transmit/receive switch 154 is controlled by a signal from the pulse generator 

45 module 121 to electrically connect the RF amplifier 151 to the coil 152 during the transmit mode and to connect the 
preamplifier 1 53 during the receive mode. The transmit/receive switch 1 54 also enables a separate RF coil (for example, 
a head coil or surface coil) to be used in either the transmit or receive mode. 

In addition to supporting the polarizing magnet 140 and the gradient coils 139 and RF coil 152, the main magnet 
assembly 141 also supports a set of shim coils 156 associated with the main magnet 140 and used to correct inhomo- 

50 geneities in the polarizing magnet field. The main power supply 157 is utilized to bring the polarizing field produced by 
the superconductive main magnet 140 to the proper operating strength and is then removed. 

The NMR signals picked up by the RF coil 152 are digitized by the transceiver module 150 and transferred to a 
memory module 160 which is also part of the system control 122. When the scan is completed and an entire array of 
data has been acquired in the memory module 160, an array processor 161 operates to Fourier transform the data into 

55 an array of image data. This image data is conveyed through the serial link 1 1 5 to the computer system 1 07 where it is 
stored in the disk memory 1 1 1 . In response to commands received from the operator console 1 00, this image data may 
be archived on the tape drive 1 12, or it may be further processed by the image processor 106 as will be described in 
more detail below and conveyed to the operator console 100 and presented on the video display 118. 
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Referring particularly to Figs. 1 and 2, the transceiver 150 includes components which produce the RF excitation 
field B 1 through power amplifier 151 at a coil 152A and components which receive the resulting NMR signal induced in 
a coil 152B. As indicated above, the coils 152A and B may be separate as shown in Fig. 2, or they may be a single 
wholebody coil as shown in Fig. 1. The base, or carrier, frequency of the RF excitation field is produced under control 

5 of a frequency synthesizer 200 which receives a set of digital signals (CF) through the backplane 118 from the CPU 
module 119 and pulse generator module 121 . These digital signals indicate the frequency and phase of the RF carrier 
signal which is produced at an output 201 . The commanded RF carrier is applied to a modulator and up converter 202 
where its amplitude is modulated in response to a signal R(t) also received through the backplane 1 18 from the pulse 
generator module 121 . The signal R(t) defines the envelope, and therefore the bandwidth, of the RF excitation pulse to 

w be produced. It is produced in the module 121 by sequentially reading out a series of stored digital values that represent 
the desired envelope. These stored digital values may, in turn, be changed from the operator console 100 to enable any 
desired RF pulse envelope to be produced. The modulator and up converter 202 produces an RF pulse at the desired 
Larmor frequency at an output 205. 

The magnitude of the RF excitation pulse output through line 205 is attenuated by an exciter attenuator circuit 206 

is which receives a digital command, TA, from the backplane 118. The attenuated RF excitation pulses are applied to the 
power amplifier 151 that drives the RF coil 152A. For a more detailed description of this portion of the transceiver 122, 
reference is made to U.S. Patent No. 4,952,877 which is incorporated herein by reference. 

Referring still to Fig. 1 and 2 the NMR signal produced by the subject is picked up by the receiver coil 152B and 
applied through the preamplifier 1 53 to the input of a receiver attenuator 207. The receiver attenuator 207 further amplifies 

20 the NMR signal and this is attenuated by an amount determined by a digital attenuation signal (RA) received from the 
backplane 118. The receive attenuator 207 is also turned on and off by a signal from the pulse generator module 121 
such that it is not overloaded during RF excitation. 

The received NMR signal is at or around the Larmor frequency, which in the preferred embodiment is around 63.86 
MHz for 1 .5 Tesla and 21 .28 MHz for 0.5 Tesla. This high frequency signal is down converted in a two step process by 

25 a down converter 208 which first mixes the NMR signal with the carrier signal on line 201 and then mixes the resulting 
difference signal with the 2.5 MHz reference signal on line 204. The resulting down converted NMR signal on line 212 
has a maximum bandwidth of 125 kHz and it is centered at a frequency of 187.5 kHz. The down converted NMR signal 
is applied to the input of an analog-to-digital (A/D) converter 209 which samples and digitizes the analog signal at a rate 
of 250 kHz. The output of the A/D converter 209 is applied to a digital detector and signal processor 210 which produce 

30 16-bit in-phase (I) values and 16-bit quadrature (Q) values corresponding to the received digital signal. The resulting 
stream of digitized I and Q values of the received NMR signal is output through backplane 1 18 to the memory module 
160 where they are employed to reconstruct an image. 

To preserve the phase information contained in the received NMR signal, both the modulator and up converter 202 
in the exciter section and the down converter 208 in the receiver section are operated with common signals. More 

35 particularly, the carrier signal at the output 201 of the frequency synthesizer 200 and the 2.5 MHz reference signal at 
the output 204 of the reference frequency generator 203 are employed in both frequency conversion processes. Phase 
consistency is thus maintained and phase changes in the detected NMR signal accurately indicate phase changes 
produced by the excited spins. The 2.5 MHz reference signal as well as 5, 1 0 and 60 MHz reference signals are produced 
by the reference frequency generator 203 from a common 20 MHz master clock signal. The latter three reference signals 

40 are employed by the frequency synthesizer 200 to produce the carrier signal on output 201 . For a more detailed descrip- 
tion of the receiver, reference is made to U.S. Patent No. 4,992,736 which is incorporated herein by reference. 

Referring particularly to Fig. 3, a preferred embodiment of a pulse sequence which may be used to acquire NMR 
data according to the present invention is shown. Actually, two pulse sequences are shown, with the first being used to 
acquire NMR data for a synchronous spin motion image (A), and the second being used to acquire NMR data for a 

45 reference image (R). In the preferred embodiment, these two pulse sequences are alternated throughout the scan such 
that the corresponding views (i.e. phase encodings) in the data sets (A) and (R) are acquired at substantially the same 
moment in time. 

The pulse sequences are fundamentally a 2DFT pulse sequence using a gradient recalled echo. Transverse mag- 
netization is produced by a selective 90° rf excitation pulse 300 which is produced in the presence of a slice select 

so gradient (G z ) pulse 301 and followed by a rephasing gradient pulse 302. A phase encoding gradient (G y ) pulse 304 is 
then applied at an amplitude and polarity determined by the view number of the acquisition. A read gradient (G x ) is 
applied as a negative dephasing lobe 306, followed by a positive readout gradient pulse 307. An NMR echo signal 309 
is acquired 40 msecs. after the rf excitation pulse 300 during the readout pulse 307 to frequency encode the 256 digitized 
samples. The pulse sequence is concluded with spoiler gradient pulses 312 and 313 along read and slice select axes, 

55 and a rephasing gradient pulse 311 is applied along the phase encoding axis (G y ). As is well known in the art, this 
rephasing pulse 31 1 has the same size and shape, but opposite polarity of the phase encoding pulse 304. The pair of 
pulse sequences are repeated 128 times with the phase encoding pulse 304 stepped through its successive values to 
acquire a 128 by 256 array of complex NMR signal samples that comprise the data set (A) and a 128 by 256 array of 
complex NMR signal samples that comprise the reference data set (R). 
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To practice the present invention an alternating magnetic field gradient is applied after the transverse magnetization 
- is produced and before the NMR signal is acquired. In the preferred embodiment illustrated in Fig. 3, the read gradient 
(Gx) is used for this function and is alternated in polarity to produce five bipolar, gradient waveforms 315. The alternating 
gradient 315 has a frequency of 200 Hz and a duration of 25 msecs. At the same time, the pulse generator module 121 

5 produces sync pulses as shown at 31 7, which are also at a frequency of 200 Hz and have a specific phase relationship 
with the alternating gradient pulses 31 5. As explained above, these sync pulses 317 activate the transducer 1 30 to apply 
an oscillating stress 31 9 to the patient which has the same frequency and phase relationship. To insure that the resulting 
waves have time to propagate throughout the field of view, the sync pulses 31 7 may be turned on well before the pulse 
sequence begins, as shown in Fig. 3. 

10 The phase of the NMR signal 309 A acquired during the first pulse sequence (A) is indicative of the movement of 
the spins. If the spins are stationary, the phase of the NMR signal is not altered by the alternating gradient pulses 315, 
whereas spins moving along the read gradient axis (x) will accumulate a phase proportional to their velocity. Spins which 
move in synchronism and in phase with the alternating magnetic field gradient 215 will accumulate maximum phase of 
one polarity, and those which move in synchronism, but 180° out of phase with the alternating magnetic field gradient 

15 215 will accumulate maximum phase of the opposite polarity. The phase of the acquired NMR signal 309 A is thus affected 
by system phase errors and random movement of spins along the x-axis, as well as the "synchronous" movement of 
spins along the x-axis. 

The reference pulse sequence is designed to measure the signal phase produced by sources other than synchro- 
nized spin movement. This is accomplished by repeating the identical pulse sequence, but without applying the oscillating 

20 stress 319. As a result, the phase of the acquired NMR signal 309 R will be affected by "static" system phase errors 
caused by field inhomogeneities and the like as well as the phase due to random spin movement along the x-axis. 
However, there will not be a phase component due to synchronous spin movement and the reference phase <|> R can, 
therefore, be subtracted from the phase (|) A to yield the phase ((|)) due solely to synchronous spin motion. 

The pulse sequence in Fig. 3 can be modified to measure synchronous spin movement along the other gradient 

25 axes (y and z). For example, the alternating magnetic field gradient pulses may be applied along the phase encoding 
axis (y) as indicated by dashed lines 321 , or they may be applied along the slice select axis (z) as indicated by dashed 
lines 322. Indeed, they may be applied simultaneously to two or three of the gradient field directions to "read" synchronous 
spin movements along any desired direction. 

The present invention may be implemented using most types of MR imaging pulse sequences. Gradient echo 

30 sequences can be readily modified to incorporate the alternating gradient as illustrated in the preferred embodiment. In 
some cases, however, the characteristics of a gradient echo sequence may not be ideal for a particular application of 
the technique. For example, some tissues (such as those with many interfaces between materials with dissimilar mag- 
netic susceptibilities) may have a relatively short T2* relaxation time and therefore may not provide enough signal to 
obtain a noise-free image at the required echo delay time. In this setting, a spin echo implementation of the invention 

35 may be ideal, because for a given echo delay time TE, this pulse sequence is much less sensitive to susceptibility effects 
than a gradient echo sequence. When a spin echo pulse sequence is used, the alternating magnetic field gradient can 
be applied either before and/or after the 1 80° rf inversion pulse. However, if the alternating gradient is applied both before 
and after the rf inversion pulse, the phase of the alternating magnetic field gradient must be inverted 180° after the rf 
inversion pulse in order to properly accumulate phase. 

40 In other applications, reduced acquisition time may be desirable. Fast spin echo and RARE sequences are rapid 
MRI sequences that acquire multiple views per TR cycle by applying different phase encoding gradients to each echo 
in a spin echo train. If 1 6 echoes, for instance, are acquired in each repetition of the sequence, then the total acquisition 
time for a complete image will be reduced by a factor of 16. One approach for modifying a fast spin echo sequence to 
implement the present invention is to insert the alternating gradient pulses between the initial 90° RF pulse and the first 

45 180° RF refocussing pulse, followed by a similar but inverted set of gradient pulses. The first echo in the train might be 
at a TE of 40-60 msec, but the spacing between subsequent echoes could be as short as 12-15 msec. 

Echo-planar imaging ("EPI") is another approach for high speed MR acquisition. In one version of this technique, 
the spin echo created by standard 90° and 180° RF pulses is broken up into a series of 64-128 short gradient echoes 
by rapidly reversing the readout gradient. A different phase encoding is applied to each of the gradient echoes and 

50 therefore the acquired data from one shot of the sequence can, in principle, be used to reconstruct a complete image. 
One approach for modifying such an echo-planar sequence to implement the present invention is to insert alternating 
gradient pulses between the initial 90° RF pulse and the 180° RF refocussing pulse, followed by a similar but inverted 
set of alternating gradient pulses. Such a sequence permits phase images to be obtained in only a few seconds or less. 
The number of cycles of the alternating magnetic field gradient used in each pulse sequence depends on the strength 

55 of the applied gradient field, the frequency of the synchronous movement to be measured, and the TE time of the pulse 
sequence. The phase sensitivity of the pulse sequence to synchronous spin movement is proportional to the integrated 
product of alternating gradient field amplitude and the displacement over time. The sensitivity may be increased by 
increasing the amplitude of the gradient field pulses and by increasing the area under each pulse by making them as 
"square" as possible. The duration of each gradient pulse is limited by the desired synchronous frequency, and hence 
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more cycles of the alternating gradient waveform are required at higher frequencies to produce the same sensitivity as 
a lower frequency alternating gradient of the same amplitude and wave shape. 

Phase sensitivity to synchronous motion can also be increased by applying both the alternating gradient pulses 315 
and the sync pulses 31 7 during the reference pulse sequence. However, when this is done the phase of the alternating 

5 magnetic field gradient 315 must be inverted 180° relative to the sync pulses 317 so that the sign of the accumulated 
phase is reversed. In addition, the magnetic field gradients should be flow compensated as described, for example, in 
U.S. patent No. 4,728,890 by Pattany et al., which is hereby incorporated by reference. Such flow compensation removes 
any phase component due to random spin motion without significantly affecting the sensitivity to synchronous spin 
motion. Consequently, when the phase difference image is produced the phase accumulations due to synchronous spin 

10 motion add together, while phase accumulations due to other sources subtract and are thereby nulled. 

If the synchronous frequency to be measured exceeds the frequency at which the magnetic field gradient can be 
switched, a lower harmonic frequency may be used for the alternating gradient. As long as the time period of one cycle 
of the alternating gradient corresponds with an odd number of cycles of the synchronous spin motion, the phase of 
synchronous spin motion will accumulate. For example, the alternating gradient frequency may be 1/3, 1/5, 1/7, 1/9, etc. 

is of the synchronous motion frequency. 

The oscillating stress may be applied by the transducer 130 in a number of ways. By starting the sync pulses 317 
well before the alternating magnetic field gradient 315 as shown in Fig. 3, the synchronous spin motion propagates 
throughout the field of view of the reconstructed image. This will image the steady-state conditions in the medium when 
the oscillating stress is applied. If the sync pulses 31 7 are turned off just before the alternating gradient 31 5 is applied, 

20 spins adjacent to the transducer 130 are moving with less amplitude or not at all during the phase accumulation time 
period. This may be desired, for example, when regions deep beneath the surface are of primary interest and large 
strain effects in the image near the transducer 130 can be suppressed. This is shown, for example, in Fig. 10a where 
the synchronous stress variations are depicted by the waveforms 336 and 338. The stress oscillations may be terminated 
one cycle short, as shown by waveform 336 when phase accumulation near the transducer 130 is to be reduced only 

25 slightly, or it may be terminated sooner as shown by waveform 338 if more phase accumulation is to be prevented. 

Referring particularly to Fig. 10b, the synchronous stress variations may only be applied for a moment in order to 
observe the mechanical impulse response or transient response of the gyromagnetic material. As indicated by waveforms 
340-343, the burst of stress oscillations may be applied at different times (to, ti, t 2 , t3, etc.) before the application of the 
alternating magnetic field gradient pulses 315 and a two-dimensional strain image reconstructed at each timing. A three- 

30 dimensional data set is thus acquired in which one of the dimensions is the time at which the stress oscillations are 
applied. The response of the spins at any image pixel location in the field of view can then be plotted as a function of 
time. The same procedure may be used to measure the response of the medium when the oscillating stress is turned 
off, as shown by the waveforms 344-346 in Fig. 10b. 

In the preferred embodiment the frequency of the sync pulses 317 and the alternating magnetic field gradient 315 

35 are constant during each scan. Other information can be acquired, however, if the frequency of either is varied. For 
example, if the frequency of the alternating gradient 315 is swept through a series of frequencies during each pulse 
sequence, or a series of acquisitions of each view are acquired at different alternating gradient frequencies, the mechan- 
ical response of the gyromagnetic material to a band of stress frequencies can be measured. Similarly, the sync pulses 
317 can be swept through a series of frequencies during each pulse sequence or during a series of acquisitions of the 

40 same view to produce a band of synchronous spin movement frequencies. The same result may also be achieved by 
changing the shape of the applied oscillatory stress from substantially sinusoidal as depicted in Fig. 3 to some other 
shape such as a square wave, which contains higher frequency harmonics. 

Such multi-frequency excitations and detections may be required in some applications where reflections cause 
standing waves. One such application, for example, is imaging the head where standing waves are formed by reflections 

45 off the skull. By imaging with multiple frequencies and adding or averaging the resulting signals, the effects of such 
standing waves can be reduced. 

A scan using the pulse sequence of Fig. 3 is carried out under the direction of a program executed by the NMR 
system of Fig. 1 . Since many different mechanical properties can be imaged according to the present invention, many 
different programs may be stored and called up for use by the operator. Some of these programs are illustrated in the 

so drawings by flow charts and described in detail below. However, it can be appreciated that many variations from these 
specific examples are possible. 

Referring particularly to Fig. 4, a scan may be performed according to the present invention to acquire NMR data 
from which ten images are reconstructed showing the propagation of strain through the gyromagnetic medium. The 
program for this scan is entered at 400 and the pulse sequence of Fig. 3 is downloaded to the pulse generator module 

55 121 . The sync pulses 21 7 in this pulse sequence are timed to be in phase with the alternating motion encoding gradient 
21 5 as indicated at process block 402. The pulse sequence is then performed the necessary number of times to acquire 
two complete NMR data sets as indicated at process block 404. As described above, one of these data sets (A) is phase 
encoded with the transducer 1 30 energized, and the other data set (R) is a reference that is acquired with the transducer 
130 deenergized. As will be described in detail below with reference to Figs. 6a-6c, a single phase image (<|)) is then 
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reconstructed from these two NM R data sets (A) and (R) as indicated at process block 406. The intensity of each pixel 
in this phase image W) is determined by the phase imparted to spins moving in synchronism with the applied oscillatory 
stress. This synchronous movement is a measure of the strain in the gyromagnetic material and the phase image ((|)) 
is, therefore, an image which shows by the brightness of its pixels the strain at corresponding locations (i.e. voxels) in 
5 the image plane. 

This phase image is a snapshot of the strain at one moment in the applied stress cycle, and to see the strain at 
other moments, the phase relationship of the sync pulses 217 is changed, as indicated at process block 408, and the 
acquisition and image reconstruction process is repeated. The scan continues to acquire NMR data and reconstruct 
strain images at successive increments (for example, 36°) of sync pulse phase. When an entire 360° cycle has been 
10 captured in successive strain images, the system branches at decision block 410. The set of strain images may then be 
viewed consecutively to "animate" the propagation of strain through the gyromagnetic medium as indicated at process 
block 412. 

The single or "animated" strain images obtained by this technique provide a view of the propagation of strain waves 
through the tissues imaged. Even without further processing, it is possible to qualitatively characterize tissue and to 

is detect abnormalities by noting various features of the depicted waves such as wave-spacing, spatial phase, diffraction, 
and attenuation. The animated images are particularly useful for detecting small lesions such as tumors that differ in 
compliance from surrounding tissues, because the localized phase shift is more easily identified than in a single static 
strain image. A simple approach for integrating the information provided by these strain images with conventional mod- 
ulus (anatomy) images provided by the same pulse sequence, is to multiply the corresponding gray scale values of the 

20 two types of images, thereby obtaining a composite image which provides both types of information. 

Although the primary focus of most medical applications of the invented imaging technique is to characterize tissue 
and to detect focal abnormalities such as tumors, the technique is also useful for characterizing acoustic systems per 
se. For instance, the technique may be used to characterize an array of acoustic energy sources that are intended for 
tissue ablation or lithotripsy. By exciting each acoustic driver individually and obtaining a strain image according to the 

25 present invention, the precise strain pattern associated with each transducer can be shown in situ, thereby showing the 
effects of diffraction and phase shifts caused by intervening tissue. These strain images could then be used to determine 
the proper phasing pattern to be applied to each of the elements in the acoustic drive array, so as to yield maximum 
strain at the target within the patient. As a final check, the tailored excitation so-produced could be applied to all of the 
drivers in the array at a low power level and a strain image could be obtained to confirm that the acoustic strain is 

30 maximized at the target tissue. The strain image could also be used to exclude the possibility of collateral damage to 
other locations when therapeutic power levels are applied to the acoustic drivers. 

While it is apparent that there are many medical applications of the invented imaging technique, it is also clear that 
the strain images provided by the technique are useful for a variety of industrial design and testing applications. For 
instance, it is a common requirement to assess the vibrational characteristics and stability of complex engine parts. 

35 Because it has not been possible to readily measure the internal strain caused by vibrational loading on such parts, it 
is often necessary to simulate such parameters using finite element analysis. The invented imaging method could be 
used to directly observe these phenomena. In this application, a model of the part would be made from a porous material 
that has a scaled approximation of mechanical properties of the real part. The pores of the model part would then be 
impregnated with a material such as oil, or water that would provide an MRI signal. Vibrational loading could then be 

40 applied, and the internal pattern of strain could be imaged using the invented method. As in medical applications, the 
objectives of non-medical applications include (1) characterization of material properties, (2) detection of focal flaws or 
abnormalities, (3) depicting the pattern of propagating waves in a transmitting medium, (4) characterizing an acoustical 
driver or array of drivers. 

While strain images provide valuable information concerning the mechanical properties of the gyromagnetic media, 
45 the present invention enables one to measure and image many other mechanical properties as well. A scan which 
acquires the data necessary for such measurements may be carried out under the direction of a program illustrated by 
the flow chart in Fig. 5. 

Referring particularly to Fig. 5, the scan begins at block 420 with the stress set to be applied along one axis. The 
direction of the alternating gradient is then set along the same axis at process block 422 and the phase of the sync 
so pulses relative to the motion encoding gradient is set to zero as indicated at process block 424. The pulse sequence of 
Fig. 3 is then executed to acquire the two image data sets (A) and (R) at process block 426, as described above. A 
phase difference image (<|)) is then produced as indicated at process block 428, and the system then loops back at 
decision block 430. 

During the second pass, the sync pulse phase is offset 90° from the phase of the alternating gradient as indicated 
55 at process block 432, and another phase image (<|)) is acquired and reconstructed. As a result, after two passes as 
detected at decision block 430, two phase images indicative of strain S L and S Q are produced as defined above in 
respective equations (7) and (9). As discussed above, considerable information about the mechanical properties of the 
gyromagnetic media can be calculated from these two strain images alone. 
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However, the scan may proceed as indicated at decision block 432 to acquire further data from which more exhaus- 
tive mechanical properties can be measured. More specifically, the direction of the alternating gradient is changed to 
align along an orthogonal axis as indicated at process block 434, and two more strain images are acquired and recon- 
structed. This may be repeated again as determined at decision block 432 with the alternating gradient switched in 

5 direction to a third orthogonal axis at process block 434. Thus, two strain images S L and S Q are produced with the 
alternating gradient directed along each of the three coordinate axes x, y and z to completely define the strain produced 
by the stress applied along one axis. 

As indicated at decision block 436, the scan may continue by altering the direction of the applied stress. In most 
applications, this is accomplished by switching the sync pulses to a different stress transducer 1 30 that is oriented in an 

10 orthogonal direction to the first transducer 130, although it is also possible to pause the scan at this point and instruct 
the operator to change the location of the transducer 130. In any case, the direction of the applied oscillatory stress is 
switched as indicated at process block 438 and the system loops back to process block 422 to repeat the above-described 
measurements. This may be repeated again for a third orthogonal stress axis, as determined at decision block 436. If 
the scan is exhaustively executed, the following eighteen images are acquired and reconstructed: 
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It should be apparent that not all of these images need be acquired during a scan. Instead, the operator selects the 

elastic property or properties to be imaged and the program is configured to acquire and reconstruct the necessary 
so strain images. From these strain images the desired elastic property can be calculated at each pixel, as indicated at 

process block 440 and described above under the heading "General Description of The Invention". The calculated values 

are then employed to modulate the intensity of its pixel location in an image. 

In the scans described above a phase image is produced from two acquired NMR data sets (A) and (R). Referring 

particularly to Figs. 6A-C, there are three methods which may be used for performing this function. As shown in Fig. 6A 
55 at process block 450, in the first method a two-dimensional Fourier transformation ("2DFT") is performed on each 

acquired data set (A) and (R). The phase ty' A is then calculated for each pixel location in the (A) image, and the phase 

ty R is calculated at each pixel of the (R) image as indicated at process block 452. 

rain a A /r A 
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f R = Tan 1 Qrfr R 

The phase difference {§) at each image pixel is then calculated at process block 454 by subtracting the reference phase 
ty' R from the synchronous motion encoded phase ty' A . 

5 The second method is referred to in the art as the complex difference technique. Referring particularly to Fig. 6B, 

the first step as indicated at process block 456 is to calculate the complex difference of the I and Q values of the corre- 
sponding sampled NM R signals in the two data sets (A) and (R). After these values are linearized, as indicated at process 
block 457, a two-dimensional Fourier transformation is performed on the resulting complex difference, as indicated at 
process block 458, and the phase at each pixel in the resulting image is calculated as indicated at process block 460. 

10 The third method shown in Fig. 6C is very similar to the method of Fig. 6A. As indicated at process block 462, a 
two-dimensional Fourier transformation is performed on each NMR data set (A) and (R) and the complex difference 
between the resulting two images is calculated at process block 464. The resulting values are linearized as indicated at 
process block 465 and the phase of the resulting complex difference image is calculated at each pixel as indicated at 
process block 466. 

15 It should be apparent to those skilled in the art that many variations are possible in the particular method used to 
acquire the NMR data and reconstruct an image indicative of the spin phase accumulated as a result of synchronous 
spin motion. 

In the preferred embodiments described above, images are reconstructed in which pixel intensity is determined 
directly by the phase of the synchronous motion encoded NMR signals. Other, less direct methods for detecting the 
20 synchronous motion are also possible. For example, if the frequency of the synchronous motion and alternating magnetic 
field gradient is increased, the spatial wavelength of the strain waves in normal tissues approaches the dimension of an 
imaged voxel. As a result, significant NMR signal phase dispersion occurs within each voxel, and in a conventional image 
in which the brightness of each pixel is determined by the "modulus" (i.e. 

y/l 2 * O 2 ) 



of the signal at each voxel, normal tissues will appear dark. On the other hand, the same strain waves passing through 

30 tumors and other less compliant tissues have a longer spatial wavelength and correspondingly less intravoxel phase 
dispersion. As a result, the signal intensity, or modulus, at these lesions will be much greater than the surrounding normal 
tissue and much brighter in the reconstructed image. 

Another less direct method for detecting the synchronous spin motion is to encode the resulting phase accumulation 
in longitudinal magnetization by tipping the phase encoded transverse magnetization back to the longitudinal axis. A 

35 pulse sequence which uses this synchronous phase encoded longitudinal magnetization will produce a modulus image 
that depicts the strain wave pattern in the subject. 

As indicated above, synchronous spin motion is produced using a transducer 130 that applies an oscillatory force 
to the gyromagnetic media. This force produces a corresponding stress in the gyromagnetic media which is propagated 
therethrough as strain in accordance with its elastic properties. This strain is the synchronous spin motion detected and 

40 imaged according to the present invention. 

Referring particularly to Figs. 7A and 7B, the preferred embodiment employs a transducer 1 30 suitable for applying 
an oscillatory force externally to a patient. It includes a coil of wire 510 wound on a bobbin 512, with the coil axis 514 
directed perpendicular to the polarizing magnetic field B 0 . The coil 510 is 400 turns of 30 AWG. copper wire, and its 
leads 516 are connected directly to the amplifier 129 as described above. The bobbin 512 is mounted to a flexible arm 

45 518 that is attached to, and cantilevered from a supporting block 520. When a current passes through the coil 510, the 
magnetic field which it produces interacts with the polarizing field B 0 . The bobbin 51 2 is thus twisted to bend the flexible 
arm 518 either upward or downward depending on the direction of current flow. By alternating the direction of current 
flow, the coil 510 twists back and forth to produce a corresponding alternating force which causes the flexible arm 518 
to oscillate as indicated by the arrow 522. A strap 521 extends through an opening in the supporting block 520 and 

50 securely fastens it to the patient. 

The oscillatory motion of the arm 518 is coupled to the subject 530 by an applicator 524. In its simplest form the 
applicator 524 is comprised of a tube 526 of the desired length which is secured at one of its ends to the arm 518. A 
pressure plate 528 is secured to the other end of the tube 526. The pressure plate 528 rests on the subject 530 to be 
imaged and its oscillatory movement produces a corresponding oscillatory compressive force that generates the desired 

55 synchronous spin motion. 

Other applicators 524 may be used with the transducer 130 to couple the oscillatory motion of the flexible arm 518 
to the subject. For example, the applicator 524 may take the form of a probe which is inserted into an opening in the 
subject to rest against a particular structure (such as the prostate gland), or it may be inserted through tissue and its 
end anchored to an internal structure (such as the liver). Many shapes and sizes are possible. 
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With the exception of the coil 51 0, all elements of the transducer 1 30 are constructed using non-magnetic materials. 
Thermoplastics such as acrylic may be used, for example, because such materials can be molded or machined into 
many desired shapes and they are relatively benign to the tissues which they contact. 

For some applications it may be possible to produce the alternating gradient field pulses 315 with the same trans- 

5 ducer 1 30 used to apply stress to the subject. In this case, the transducer 1 30 includes a coil that acts as a local gradient 
coil for surrounding tissues. The same alternating current applied to this coil to produce the alternating magnetic field 
gradient pulses 315 also cause the coil to move in an oscillatory manner at the same frequency and phase. Such a coil 
might be used, for example, to detect tumors in the prostate gland. 

While application of a synchronous motion is preferably accomplished by the direct application of a force to the 

10 patient, indirect synchronous motion producing mechanisms are also possible. Electrical or other non-mechanical stim- 
ulation may be applied to the patient to produce synchronous motion in the tissues to be imaged. For example, synchro- 
nous brain function stimulation may produce synchronous motion at the cellular, or sub-cellular level. Synchronous 
electrical stimulation of selected muscles will, of course, produce synchronous motion in the muscles themselves, but 
also, detectable synchronous motion may be produced in the associated nerves and in tissues surrounding the muscles. 

is Rather than applying a separate synchronous stress or synchronous stimulation, an endogenous periodic motion 
may also be used to practice the present invention. In this case the alternating motion encoding magnetic field gradient 
must be synchronized with the endogenous motion, and for best results, this requires the detection of a signal indicative 
of the frequency and phase of the endogenous motion. For example, the rapid movement of an abnormal heart valve 
may be detected with a phonocardiogram that senses the sound produced by the valve. The resulting electrical signal 

20 is used to synchronize the alternating magnetic field gradient. The resulting synchronous motion of the valve itself, as 
well as surrounding tissues and the blood flowing through the valve may be imaged. Similar results may be achieved 
with the heart using the ECG signal. 

Contrast agents are used in many imaging modalities to enhance image contrast in a diagnostically useful manner. 
Such a strategy is also possible with the present invention when the contrast agent operates to change the elastoviscous 

25 characteristics of the tissues of interest. For example, an osmotic agent which affects mechanical properties in a manner 
that may be detected using the present invention may be used. A hyper-osmolar agent will draw water out of the target 
cells and make them more compliant, whereas a hypo-osmolar agent will add water to the target cells and make them 
rigid, or less compliant. This compliance difference can be imaged with the present invention and thus the affected cells 
can be contrasted against other cells that do not respond to the contrast agent. 

30 The present invention may also be used to enhance magnetic resonance angiograms ("MRA"). MRA relies on the 
movement of blood in the direction of one or more motion encoding magnetic field gradients to differentiate between 
stationary blood vessel walls and the blood which they carry. If the blood flow is slow, as in venous blood vessels, this 
contrast mechanism is less effective. The present invention may be used to produce an MRA by applying a compressive 
force at the synchronous frequency which sends compression waves through the vasculature. Since blood is an incom- 

35 pressible liquid, the synchronous compression wave will travel quite efficiently through the vessels as an oscillatory 
motion of the blood cells at the synchronous frequency. The alternating magnetic field gradient is "tuned" to this syn- 
chronous frequency and the resulting phase image will depict the vasculature. An example of this application of the 
invention will now be described. 

Referring particularly to Fig. 1 1 , a cuff 550 is used in place of the transducer 1 30 in this embodiment of the invention. 

40 The cuff 550 receives the synchronous drive signal produced by the amplifier 129 (Fig. 1) at leads 551 and it constricts 
around the subject's ankle. The resulting synchronous compression waves in the venous structure and the arterial struc- 
ture travel up the patient's leg 553 as indicated by arrows 555. The acquisition of a phase image as described above 
using the pulse sequence of Fig. 3 with a synchronous alternating magnetic field gradient 315 directed along the length- 
wise dimension of the leg 553, may be used to produce an angiogram. To improve the SNR of the image additional 

45 phase images with the alternating magnetic field gradient directed along the other two orthogonal axes may also be 
acquired and summed together. 

Claims 

so 1. A method for producing an image with an NMR imaging system which indicates a mechanical characteristic of a 
subject within the field of view of the NMR imaging system, the steps comprising: 

a) imparting mechanical wave motion to spins located in the subject (317, 129, 130); 

b) conducting a scan of the subject (404, 426) with the NMR imaging system using a pulse sequence (Fig. 3) 
55 in which an alternating magnetic field gradient (315) is employed to detect the wave motion of the spins; and 

c) reconstructing a phase image from NMR data (Fig. 6A-C) acquired during the scan (406, 428) and employing 
the phase information therein to indicate the mechanical characteristic of the subject at each pixel location of 
the phase image. 
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2. The method as recited in claim 1 in which the reconstructed image (406, 428) indicates the phase of the acquired 
NMR data caused by the mechanical wave motion of the spins. 

3. The method as recited in claim 1 in which the subject (530) is a gyromagnetic medium which has elastic properties 
5 and the mechanical wave motion is imparted to the spins by applying an oscillating stress (319) to a portion of the 

subject. 

4. The method as recited in claim 3 in which the subject is a portion of a living animal. 

10 5. The method as recited in claim 1 , 2 or 3 in which the mechanical characteristic is one of a group of characteristics 
(440) such as strain, Young's modulus, Poisson's ratio , the bulk modulus, or the shear modulus. 

6. The method as recited in claim 1 in which the subject is a living animal and the mechanical wave motion is imparted 
to the spins by applying a stimulus to the living animal. 

15 

7. The method as recited in claim 2 in which the alternating magnetic field gradient (315) alternates a plurality of cycles 
before the acquisition of the NMR signal. 

8. The method as recited in claim 1 in which the mechanical wave motion is imparted to the spins by applying a stress 
20 (31 9) which varies periodically in magnitude. 

9. The method as recited in claim 2 in which a second set of NMR data is acquired (404R, 426R) from the subject 
without imparting mechanical wave motion to the spins, and the phase indicated by the second set of NMR data is 
subtracted from the phase indicated by said reconstructed image (454, 456, 464). 

25 

10. The method as recited in claim 2 in which a second phase image is produced from the subject by repeating (430) 
steps a), b) and c) with the relative phase between the mechanical wave motion and the alternating magnetic field 
gradient (315) being different than that used to produce said reconstructed image, and producing an image indicative 
of a mechanical characteristic of the subject (440) by combining said second phase image with said reconstructed 

30 image. 

11. The method as recited in claim 2 in which a second phase image is produced from the subject by repeating (432) 
steps a), b) and c) with the alternating magnetic field gradient (321 , 322) being oriented in a different direction (434) 
than that used to produce said reconstructed image, and producing an image indicative of a mechanical character- 

35 istic of the subject (440) by combining the second phase image with said reconstructed image. 

1 2. The method as recited in claim 1 1 in which the respective alternating magnetic field gradients (31 5, 321 , 332) used 
to acquire NMR signals for the two images are orthogonal. 

40 13. The method as recited in claim 10 in which the mechanical characteristic (440) is one of a group of characteristics 
such as compliance, the velocity at which a strain wave is propagated in the subject, or the Young's modulus of the 
subject. 

14. An NMR system for producing an image which indicates a mechanical characteristic of a subject placed within the 
45 field of view of the NMR system, which comprises: 

a transducer (1 30) coupled to the subject for imparting mechanical wave motion to spins located in the subject; 
a magnet (140) for producing a polarizing magnetic field in the field of view; 
an RF transmitter (150) for producing transverse magnetization in the spins; 
an RF receiver (150) for receiving NMR signals produced by the spins; 
so magnetic field gradient means (1 39) for producing a magnetic field in the field of view having a gradient in a 

first direction; 

a pulse generator (121) for producing a pulse sequence in which the RF transmitter produces transverse 
magnetization in the spins, the magnetic field gradient means produces an alternating magnetic field gradient that 
detects the wave motion of the spins, and the RF receiver receives NMR signals produced by the spins; and 
55 an image ^constructor (106) for reconstructing a phase image from the received NMR signals and which 

employs the phase information therein to indicate the mechanical characteristic of the subject at each pixel location 
of the phase image. 
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15. The NMR system as recited in claim 14 in which the transducer imparts the mechanical wave motion by applying 
a stress to the subject which varies periodically in magnitude. 

16. The NMR system as recited in claims 14 or 15 in which the subject is a gyromagnetic medium which has elastic 
properties and the transducer imparts the mechanical wave motion to the spins by applying an oscillating stress to 
a portion of the subject. 

17. The NMR system as recited in any of the claims 14 to 16 in which the image reconstructor employs the phase of 
the received NMR signals to indicate the mechanical characteristic of the subject. 

18. The NMR system as recited in any of the claims 14 to 1 7 in which the magnetic field gradient means alternates the 
polarity of the magnetic field gradient a plurality of cycles before receipt of the NMR signals. 

19. The NMR system as recited in claim 17 in which the pulse generator produces a second pulse sequence while the 
transducer is turned off to receive second NMR signals, and the image reconstructor also employs the phase of the 
received second NMR signals to indicate the mechanical characteristic of the subject. 

20. The NMR system as recited in claim 1 7 in which the pulse generator produces a second pulse sequence and second 
NMR signals are received in which the relative phase between the mechanical wave motion imparted by the trans- 
ducer and the alternating magnetic field gradient is different than the relative phase used for receiving said NMR 
signals, and the image reconstructor also employs the phase of the received second NMR signals to indicate the 
mechanical characteristic of the subject. 

21. The NMR system as recited in claim 17 in which the pulse generator produces a second pulse sequence in which 
the magnetic field gradient means produces an alternating magnetic field gradient oriented in a different direction 
than that produced during said pulse sequence, and second NMR signals are received, and in which the image 
reconstructor also employs the phase of the received second NMR signals to indicate the mechanical characteristic 
of the subject. 

22. The NMR system as recited in claim 21 in which the orientation of the alternating magnetic field gradient produced 
during the second pulse sequence is orthogonal to the orientation of the alternating magnetic field gradient produced 
during said pulse sequence. 
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Description 

Background of the Invention 

5 [0001] The field of the invention is nuclear magnetic resonance imaging methods and systems. More particularly, 
the invention relates to the enhancement of MR image contrast. 

[0002] The physician has many diagnostic tools at his or her disposal which enable detection and localization of 
diseased tissues. These include x-ray systems that measure and produce images indicative of the x-ray attenuation 
of the tissues and ultrasound systems that detect and produce images indicative of tissue echogenicity and the bound- 

10 aries between structures of differing acoustic properties. Nuclear medicine produces images indicative of those tissues 
which absorb tracers injected into the patient, as do PET scanners and SPECT scanners. And finally, magnetic reso- 
nance imaging ("MRI") systems produce images indicative of the magnetic properties of tissues. It is fortuitous that 
many diseased tissues are detected by the physical properties measured by these imaging modalities, but it should 
not be surprising that many diseases go undetected. 

15 [0003] Historically, one of the physician's most valuable diagnostic tools is palpation. By palpating the patient a 
physician can feel differences in the compliance of tissues and detect the presence of tumors and other tissue abnor- 
malities. Unfortunately, this valuable diagnostic tool is limited to those tissues and organs which the physician can feel, 
and many diseased internal organs go undiagnosed unless the disease happens to be detectable by one of the above 
imaging modalities. Tumors (e.g. of the liver) that are undetected by existing imaging modalities and cannot be reached 

20 for palpation through the patient's skin and musculature, are often detected by surgeons by direct palpation of the 
exposed organs at the time of surgery. Palpation is the most common means of detecting tumors of the prostate gland 
and the breast, but unfortunately, deeper portions of these structures are not accessible for such evaluation. An imaging 
system that extends the physician's ability to detect differences in tissue compliance throughout a patient's body would 
extend this valuable diagnostic tool. 

25 [0004] Any nucleus which possesses a magnetic moment attempts to align itself with the direction of the magnetic 
field in which it is located. In doing so, however, the nucleus precesses around this direction at a characteristic angular 
frequency (Larmor frequency) which is dependent on the strength of the magnetic field and on the properties of the 
specific nuclear species (the magnetogyric constanty of the nucleus). Nuclei which exhibit this phenomena are referred 
to herein as "spins", and materials which contain such nuclei are referred to herein as "gyromagnetic". 

30 [0005] When a substance such as human tissue is subjected to a uniform magnetic field (polarizing field B 0 ), the 
individual magnetic moments of the spins in the tissue attempt to align with this polarizing field, but precess about it 
in random order at their characteristic Larmor frequency. A net magnetic moment M z is produced in the direction of 
the polarizing field, but the randomly oriented magnetic components in the perpendicular, or transverse, plane (x-y 
plane) cancel one another. If, however, the substance, or tissue, is subjected to a magnetic field (excitation field B A ) 

35 which is in the x-y plane and which is near the Larmor frequency, the net aligned moment, Mz, may be rotated, or 
"tipped", into the x-y plane to produce a net transverse magnetic moment M t , which is rotating, or spinning, in the xy 
plane at the Larmor frequency. The practical value of this phenomenon resides in the signal which is emitted by the 
excited spins after the excitation signal B 1 is terminated. There are a wide variety of measurement sequences in which 
this nuclear magnetic resonance ("NMR") phenomena is exploited. 

40 [0006] When utilizing NMR to produce images, a technique is employed to obtain NMR signals from specific locations 
in the subject. Typically, the region which is to be imaged (region of interest) is scanned by a sequence of NMR meas- 
urement cycles which vary according to the particular localization method being used. The resulting set of received 
NMR signals are digitized and processed to reconstruct the image using one of many well known reconstruction tech- 
niques. To perform such a scan, it is, of course, necessary to elicit NMR signals from specific locations in the subject. 

45 This is accomplished by employing magnetic fields (G x , G y , and G z ) which are superimposed on the polarizing field 
B 0 , but which have a gradient along the respective x, y and z axes. By controlling the strength of these gradients during 
each NMR cycle, the spatial distribution of spin excitation can be controlled and the location of the resulting NMR 
signals can be identified. 

[0007] It is well known that NMR can be used to detect and image the movement of spins. As disclosed in U.S. patent 
50 No. Re. 32,701 entitled "NMR Scanner With Motion Zeugmatography", acquired NMR signals can be sensitized to 
detect moving spins by applying a bipolar magnetic field gradient at the proper moment in each NMR measurement 
sequence. The phase of the resulting NMR signal measures the velocity of spins along the direction of the motion 
sensitizing magnetic field gradient. With more complex motion sensitizing magnetic field gradients, higher orders of 
motion, such as acceleration and jerk can also be measured with this method. 
55 [0008] The document Radiology, vol. 191, no.3, June 1994, p. 835 to 839, J.E. Drace et al "Elastic Deformation on 
Tendons and Myotendinous Tissue:Measurement by Phase-Contrast MR imaging", discloses a method for calculating 
strain on a tendon by acquiring a series of velocity NMR images as the tendon changes shape due to being subjected 
to cyclical translation, rotation and stretch. Strain is calculated by observing differences in velocity at different points 
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in the tendon. The tendon is stretched in the same manner as each view is acquired. The motion cycle is repeated 
once for each amplitude of the spatial phase-encoding gradient. 

[0009] The method disclosed in the document "Visualizing Tissue Compliance With MR", Society of Magnetic Res- 
onance, Second Meeting, 1994, Abstract, August 1994, D.B. Plewes et al, is rather similar. A subject is compressed 
5 in the same manner as each view is acquired. A compressive force is applied to the subject and then a single pulse 
sequence is performed to measure the resulting displacement. This scheme is repeated to get all views for an image. 
The pulse sequence is synchronized with the force application in the sense that the timing is the same each repetition. 
So the subject moves the same amount as each view is acquired. 

[0010] The document "Oscillatory flow in the cochlea visualized by magnetic resonance imaging technique", Proc. 

10 Natl. Acad. Sci. USA, Vol. 90, pp. 1 595-1 598, Feb. 1 993, Biophysics, Denk et al, is directed to visualize motion of fluid 
within the cochlea of the inner ear. For this purpose, oscillating magnetic field gradients are phase-locked to an external 
stimulus to generate a magnetic resonance image of fluid movement that is correlated with the stimulus in order to 
selectively visualize and quantify oscillatory fluid motion. Two data sets are acquired, one of them taken with reversed 
stimulus phase, and the data sets are subtracted in order to eliminate all signals but those due to spins moving in 

15 synchrony with the stimulus. The images derived in this document from a complex difference data set are: a magnitude 
display of the sum of an antiphase image pair and the imaginary parts of difference images. These images allow to 
produce a map of mechanical response in the cochlea and to obtain a view of cochlear fluid mechanics. The document 
also mentions the possibility of applying the method to the study of the spatial distribution of mechanical properties of 
materials. 

20 

Summary of the Invention 

[0011] The present invention relates to the imaging of mechanical characteristics of a subject's tissue and more 
particularly, the production of an NMR image in which the intensity pattern reveals the mechanical properties of the 

25 tissue and the pattern of propagation of mechanical strain waves within it. The invention concerns a method defined 
in claim 1 and a corresponding system as defined in claim 13. The further claims define additional features. The in- 
vention includes the acquisition of NMR imaging data using an NMR pulse sequence in which a magnetic field gradient 
is alternated in polarity to phase encode moving spins, and a stimulus is applied to the object being imaged and is 
oscillated in magnitude in synchronism with the alternating magnetic field gradient. A phase image is reconstructed 

30 from the acquired NMR imaging data which depicts the synchronous movement of spins throughout the object resulting 
from the applied stimulus. Mechanical characteristics of the object are revealed by analysis of such images. 
[0012] A general object of the invention is to produce an image which reveals the mechanical properties of tissues. 
The applied stimulus may be a mechanical disturbance which sets up a traveling wave in the gyromagnetic medium 
that causes cyclic strain, or displacements of the spins. By synchronizing the alternating magnetic field gradient with 

35 the oscillatory applied stress, the movement of the spins in relation to the applied force is revealed by the phase of the 
acquired NMR signals. At locations where the spins move in phase with the applied stress the phase signal is maximum, 
at one polarity, and at locations where the spins oscillate out of phase with the applied stress the phase signals are 
maximum at the opposite polarity. These phase signals may be employed to modulate pixel intensity in a reconstructed 
image and the strain due to the applied stress throughout the imaged subject's tissues is seen as a wave pattern. 

40 Changes in this wave pattern reveal the boundaries of regions of differing mechanical properties within the subject's 
tissues. 

[0013] Another object of the invention is to produce a series of images which reveal the strain wave propagation 
through asubject's tissue. By synchronizing the application of the oscillating applied stress with the alternating magnetic 
field gradient at a series of different phase relationships, a corresponding series of phase images may be reconstructed 

45 showing the movement of the strain wave pattern through the tissue. 

[0014] Yet another object of the invention is to produce an image which indicates the strain in tissues produced by 
an applied force. A first phase image is produced with the oscillating stress applied, and a second phase image is 
produced without the oscillating stress applied. The difference between these two phase images is calculated and 
provides an image indicative of the strain in the tissue. These strain images provide a measure of a fundamental 

50 quantity of the medium in the form of displacement amplitude as a strain wave propagates through it. As a result, it 
becomes possible to calculate a number of fundamental mechanical properties pertaining to the viscoelastic nature of 
the gyromagnetic medium. They include, for example, velocity of propagation, Young's modulus and other moduli. 
[0015] Another object of the invention is to produce images indicative of viscoelastic properties of tissues. By pro- 
ducing two phase images with the oscillating stress synchronized with the alternating magnetic field gradient at two 

55 different phase relationships, the wavelength of the propagated strain can be calculated at each image pixel. From 
this, the velocity of propagation can be calculated at each image pixel as well as Young's modulus. 
[0016] Another object of the invention is to produce an image which indicates Poisson's ratio throughout a subject's 
tissue. A first phase image is produced with the alternating gradient applied in the same direction as the oscillating 
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stress, and a subsequent phase image is produced with the alternating gradient applied in an orthogonal direction. 
Poisson's ratio may be estimated at each pixel location from the resulting phase images. Since the velocity and wave- 
length of the propagated strain may also be measured, the viscous properties of the gyromagnetic medium can also 
be estimated. Using Young's modulus and Poisson's ratio, other moduli such as the shear modulus and the bulk modulus 

5 can be calculated and imaged. 

[0017] Another general object of the invention is to produce an image of a patient which indicates the compliance of 
tissues. By acquiring a plurality of phase images of the patient with different arrangements of the oscillatory stress and 
the alternating magnetic field gradient, an image indicative of tissue compliance can be produced. This image provides 
the clinician with information comparable to that obtained by manually palpating the tissues. 

w [0018] Yet another object of the invention is to produce an NMR image in which contrast is enhanced for spins having 
oscillatory motion. By synchronizing the alternating magnetic field gradient with oscillatory movement of spins in the 
object being imaged, a phase image may be reconstructed in which the signals from the oscillatory spins are enhanced 
relative to stationary spins and randomly moving spins. 

[0019] The foregoing and other objects and advantages of the invention will appear from the following description. 
15 in the description, reference is made to the accompanying drawings which form a part hereof, and in which there is 
shown by way of illustration a preferred embodiment of the invention. Such embodiment does not necessarily represent 
the full scope of the invention, however, and reference is made therefore to the claims herein for interpreting the scope 
of the invention. 

20 Brief Description Of The Drawings 

[0020] 

Fig. 1 is a block diagram of an NMR system which employs the present invention; 
25 Fig. 2 is an electrical block diagram of the transceiver which forms part of the NMR system of Fig. 1 ; 

Fig. 3 is a graphic representation of a pulse sequence performed by the NMR system of Fig. 1 to practice the 
preferred embodiment of the invention; 

Fig. 4 is a flow chart which illustrates how the NMR system of Fig. 1 reconstructs a series of strain images in 
accordance with a first preferred embodiment of the invention using NMR data acquired with the pulse sequence 
30 of Fig. 3; 

Fig. 5 is a flow chart which illustrates how the NMR system of Fig. 1 produces images indicative of the elastic 
properties in accordance with a second preferred embodiment of the invention using NMR data acquired with the 
pulse sequence of Fig. 3; 

Figs. 6a-c are flow charts of three different methods used to produce phase difference images in the methods of 
35 Figs. 4 and 5; 

Figs. 7a and 7b are elevation and top views respectively of a transducer that is used in the NMR system of Fig. 1 
to practice the present invention; 

Fig. 8 is a graphic illustration of the phase accumulation in an NMR signal acquired with the pulse sequence of Fig. 3; 
Fig. 9 is a graphic illustration of phase accumulation with the alternating gradient waveform reversed 180°; 
40 Figs. 1 0a and 1 0b are graphic illustrations of alternative methods for applying the oscillating stress relative to the 

alternating gradient waveform in the pulse sequence of Fig. 3; and 

Fig. 11 is a pictorial representation of a human leg illustrating an alternative method for stimulating synchronous 
spin motion in blood vessels. 

45 General Description Of The Invention 

[0021] In general, the physical properties of tissue are measured by applying a stress (e.g. tension, pressure, or 
shear) and observing the resulting strain (e.g. elongation, compression, rotation). By measuring the resulting strain, 
elastic properties of the tissue such as Young's modulus, Poisson's ratio, the shear modulus, and the bulk modulus, 
50 can be calculated. By applying the stress in all three dimensions and measuring the resulting strain, the elastic prop- 
erties of the tissue can be completely defined. 

[0022] By observing the rate at which the strain decreases as a function of distance from the stress producing source, 
the attenuation of the strain wave can be estimated. From this, the viscous properties of the gyromagnetic medium 
may be estimated. The dispersion characteristics of the medium can be estimated by observing the speed and atten- 
55 uation of the strain waves as a function of their frequency. Dispersion is potentially a very important parameter for 
characterizing tissues in medical imaging applications. 

[0023] The present invention provides a means for measuring the strain in gyromagnetic materials such as tissues 
using NMR methods and apparatus. The Larmor equation is given by 
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co = yB (1) 

where co is the angular resonant frequency of the NM R signal produced by spins having a characteristic gyromagnetic 
ratio of y when placed in a magnetic field having a density and direction B. The magnetic field vector B can be broken 
down into two components 

B =B 0 + rG r (2) 

where B 0 is the polarizing magnetic field, r is the location of the spins, and G r is the magnetic field gradient. 

[0024] Since the angular frequency of the NMR signals produced by the spins is the rate of change of their phase, 

the phase of the spin signals as a function of time is as follows: 

(K0 = M t)dt. (3) 

[0025] Substituting equation (1 ), the relationship between NMR signal phase and the applied gradient field is obtained 

^t) = y!G(r)r(f)dt\ (4) 

where G and r are expressed as functions of time (V) for obtaining a general expression of the Larmor equation. This 
equation indicates that the NMR signal produced by moving spins will accrue a phase shift relative to that accrued by 
static spins when in the presence of a magnetic field gradient. 

[0026] If an oscillating stress is applied to tissue along the direction r at an angular frequency co p , a wave is launched 
and spins are displaced by amounts determined by the elastic properties of the tissue. If it is assumed that this prop- 
agation occurs without loss, the displacement (A) of spins at location (r) may be expressed as follows: 

A = A 0 cos (co p t + kr +9), (5) 

where A 0 is the maximum displacement produced by the applied stress, k is the wave number, and 9 is the phase offset 
of the spin displacement relative to the applied oscillating stress. The wave number k is equal to 2n radians divided 
by the wavelength (X) of the propagated wave, and if it is assumed that the spin displacement occurs for just one cycle 
(t=0 to T) of the applied oscillating stress, then the NMR signal produced by the spin will accumulate a phase indicated 
by the following expression: 



c-r 

<K c) =y / G(t') A 0 cos(o) p t+Jcr+8)dC: / (6) 

C-0 



If the magnetic field gradient G(t') is constant during this time period, no phase signal will be accumulated. However, 
if the magnetic field gradient G(t') is synchronized with the applied stress and is switched in polarity half way through 
the time period (T), the phase of the NMR signal ((|>) at the completion of the time period will be proportional to the 
displacement of the spins. This displacement along the r direction is the strain which results from the applied stress 
along the same direction r. 

[0027] The phase accumulation of the NMR signal due to spins in synchronous motion is illustrated in Fig. 8. An 
alternating magnetic field gradient that completes one bi-polar cycle during the time period T is indicated by waveform 
10. Spins in synchronous motion and perfectly in phase with the alternating magnetic field gradient 10 is illustrated by 
the waveform 12, and the resulting accumulated phase in the NMR signal produced by these spins is indicated by the 
graph 14. On the other hand, synchronous spin motion 1 80° out of phase with the alternating gradient 1 0 is indicated 
by waveform 13 and the resulting phase accumulation is shown in graph 15. In contrast, Fig. 9 illustrates at graph 16 
the phase of NMR signals produced by stationary spins, and graph 1 8 indicates the accumulated phase of constantly 
moving spins. At the completion of one cycle of the alternating gradient 1 0, the accumulated phase of stationary spins 
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is nulled and the accumulated phase of constantly moving spins is proportional to their velocity along the direction of 
the gradient. 

[0028] Because the displacements are very small, the gradient field G(t') is usually oscillated in synchronism with 
the applied stress for several cycles before the NMR signal is acquired. This enables the accumulated phase (<|>) to 
reach a significant amount. A second NMR signal is then acquired with an identical pulse sequence, but without the 
application of the oscillatory stress. This second signal provides a reference phase (<|> R ) which may be subtracted from 
the accumulated phase <|) A to yield a value indicative of longitudinal strain (S L ). 

[0029] These NMR measurements can be made, of course, with imaging gradients applied, and a strain image 
indicative of the longitudinal strain at each pixel may be reconstructed. This strain image has pixel values as follows: 

S L (t,r) = A 0 cos(co p t + kr). (7) 

[0030] If the measurement is repeated with the oscillating stress in the same direction, but with the alternating gradient 
G(t') oriented along each of the orthogonal axes, the displacements, along these axes can be determined and the 
orthogonal strain (S T ) calculated. From this information Poisson's ratio (a) can be calculated as follows: 

o = S T /S L (8) 

[0031] Further information can be learned about the elastic properties of tissues by changing the phase relationship 
of the applied oscillatory stress and the synchronized alternating motion gradient. For example, if the phases are offset 
nl2 radians, a second strain image S 0 is produced having the following pixel values: 

S 0 (t,r) = A 0 cos(co p t + kr + 

= A 0 sin(co p t + kr) (9) 

[0032] The gradient in each of the strain images S L and S 0 may be defined as follows: 

VS L = ^i + ^j (10) 

[0033] Using these gradients, the wave number (k) of the propagated wave at each image pixel can be calculated 
as follows: 

k- J\VS L \-\VS L \MVS 0 \-\VS 0 \/Js? 0 + S? L (12) 

Knowing the frequency (f) of the applied oscillatory stress, the propagation velocity (c) can then be calculated at each 
image pixel as long as the viscosity effects are not significant. 

X = 27i/k 

c = fX (13) 



If the density (p) of the gyromagnetic medium is known, the propagation velocity (c) can be used to calculate Young's 
modulus (Y): 
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Y = c 2 p. 



(14) 



[0034] Since compliance is the inverse of Young's modulus (i.e.1/Y), the compliance of the gyromagnetic medium 
may also be calculated. An image in which pixel intensity is determined by the calculated compliance has diagnostic 
value in medicine because such an image displays what a physician feels when manually palpating tissue. With the 
knowledge of Poisson's ratio (a) and the Young's modulus (Y) all other values of moduli, namely the shear modulus 
and the bulk modulus ((3) may be calculated, since only two of these four elastic properties are actually independent. 
Bulk modulus (P) can be expressed from Young's modulus and Poissons ratio, in the following way: 



Y = 3p(1 -2a). 



(15) 



Shear Modulus (ju) can be expressed as follows: 



Y = 2n(1 +a). 



(16) 



It should be apparent that the oscillatory stress may be applied in three orthogonal directions and the synchronized 
gradient field may also be applied in three separate orthogonal directions for each applied stress direction. Spin dis- 
placements (A, r|, E) in all three directions may be measured, and all the components of the strain dyadic (second- 
rank tensor field), given in the following matrix may be measured: 



dx dx dx 

dy dy dy 

3A dr\ 3E 

^ dz dz dz 



(17) 



With these measurements it is possible to characterize all the elastic properties of the gyromagnetic medium under 
investigation using the above-described calculations. 



Description Of The Preferred Embodiment 

[0035] Referring first to Fig. 1 , there is shown the major components of a preferred NMR system which incorporates 
the present invention and which is sold by the General Electric Company under the trademark "SIGNA". The operation 
of the system is controlled from an operator console 1 00 which includes a console processor 1 01 that scans a keyboard 
102 and receives inputs from a human operator through a control panel 103 and a plasma display/touch screen 104. 
The console processor 101 communicates through a communications link 116 with an applications interface module 
117 in a separate computer system 107. Through the keyboard 102 and controls 103, an operator controls the pro- 
duction and display of images by an image processor 106 in the computer system 107, which connects directly to a 
video display 1 1 8 on the console 1 00 through a video cable 1 05. 

[0036] The computer system 107 includes a number of modules which communicate with each other through a 
backplane. In addition to the application interface 117 and the image processor 106, these include a CPU module 108 
that controls the backplane, and an SCSI interface module 1 09 that connects the computer system 1 07 through a bus 
1 1 0 to a set of peripheral devices, including disk storage 1 1 1 and tape drive 1 1 2. The computer system 1 07 also includes 
a memory module 113, known in the art as a frame buffer for storing image data arrays, and a serial interface module 
114 that links the computer system 107 through a high speed serial link 115 to a system interface module 120 located 
in a separate system control cabinet 122. 

[0037] The system control 1 22 includes a series of modules which are connected together by a common backplane 
1 1 8. The backplane 1 1 8 is comprised of a number of bus structures, including a bus structure which is controlled by a 
CPU module 119. The serial interface module 120 connects this backplane 118 to the high speed serial link 115, and 
pulse generator module 121 connects the backplane 118 to the operator console 100 through a serial link 125. It is 
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through this link 125 that the system control 122 receives commands from the operator which indicate the scan se- 
quence that is to be performed. 

[0038] The pulse generator module 121 operates the system components to carry out the desired scan sequence. 
It produces data which indicates the timing, strength and shape of the RF pulses which are to be produced, and the 
5 timing of and length of the data acquisition window. The pulse generator module 121 also connects through serial link 
1 26 to a set of gradient amplifiers 1 27, and it conveys data thereto which indicates the timing and shape of the gradient 
pulses that are to be produced during the scan. 

[0039] In the preferred embodiment of the invention the pulse generator module 121 also produces sync pulses 
through a serial link 128 to a wave generator and amplifier 129. The wave generator produces a sinusoidal voltage 
10 which is synchronized to the frequency and phase of the received sync pulses and this waveform is output through a 
50 watt, dc coupled audio amplifier. A frequency in the range of 20 Hz to 1000 Hz is produced depending on the 
particular object being imaged, and it is applied to a transducer 130. The transducer 130 will be described in more 
detail below, and its structure depends on the particular anatomy being measured and imaged. In general, however, 
the transducer 1 30 produces a force, or pressure, which oscillates in phase with the sync pulses produced by the pulse 
15 generator module 121 and creates an oscillating stress in the gyromagnetic media (i.e. tissues) to which it is applied. 
[0040] And finally, the pulse generator module 121 connects through a serial link 132 to scan room interface circuit 
1 33 which receives signals at inputs 1 35 from various sensors associated with the position and condition of the patient 
and the magnet system. It is also through the scan room interface circuit 133 that a patient positioning system 134 
receives commands which move the patient cradle and transport the patient to the desired position for the scan. 
20 [0041] The gradient waveforms produced by the pulse generator module 121 are applied to a gradient amplifier 
system 127 comprised of G x , G y and G z amplifiers 136, 137 and 138, respectively. Each amplifier 136, 137 and 138 
is utilized to excite a corresponding gradient coil in an assembly generally designated 1 39. The gradient coil assembly 
1 39 forms part of a magnet assembly 141 which includes a polarizing magnet 1 40 that produces either a 0.5 or a 1 .5 
Tesla polarizing field that extends horizontally through a bore 142. The gradient coils 139 encircle the bore 142, and 
25 when energized, they generate magnetic fields in the same direction as the main polarizing magnetic field, but with 
gradients G x , G y and G z directed in the orthogonal x-, y- and z-axis directions of a Cartesian coordinate system. That 
is, if the magnetic field generated by the main magnet 140 is directed in the z direction and is termed B 0 , and the total 
magnetic field in the z direction is referred to as B z , then G x =dB z /dx, G y =dB z /dy and G z =dB z /dz, and the magnetic field 
at any point (x,y,z) in the bore of the magnet assembly 141 is given by B(x,y,z)=B 0 +G x x+G y y+G z z. The gradient mag- 
so netic fields are utilized to encode spatial information into the NMR signals emanating from the patient being scanned, 
and as will be described in detail below, they are employed to measure the microscopic movement of spins caused by 
the pressure produced by the transducer 130. 

[0042] Located within the bore 1 42 is a circular cylindrical whole-body RF coil 1 52. This coil 1 52 produces a circularly 
polarized RF field in response to RF pulses provided by a transceiver module 150 in the system control cabinet 122. 

35 These pulses are amplified by an RF amplifier 151 and coupled to the RF coil 152 by a transmit/receive switch 154 
which forms an integral part of the RF coil assembly. Waveforms and control signals are provided by the pulse generator 
module 121 and utilized by the transceiver module 150 for RF carrier modulation and mode control. The resulting NMR 
signals radiated by the excited nuclei in the patient may be sensed by the same RF coil 152 and coupled through the 
transmit/receive switch 154 to a preamplifier 153. The amplified NMR signals are demodulated, filtered, and digitized 

40 in the receiver section of the transceiver 150. The transmit/receive switch 154 is controlled by a signal from the pulse 
generator module 121 to electrically connect the RF amplifier 151 to the coil 152 during the transmit mode and to 
connect the preamplifier 153 during the receive mode. The transmit/receive switch 154 also enables a separate RF 
coil (for example, a head coil or surface coil) to be used in either the transmit or receive mode. 
[0043] In addition to supporting the polarizing magnet 140 and the gradient coils 139 and RF coil 152, the main 

45 magnet assembly 141 also supports a set of shim coils 1 56 associated with the main magnet 1 40 and used to correct 
inhomogeneities in the polarizing magnet field. The main power supply 157 is utilized to bring the polarizing field pro- 
duced by the superconductive main magnet 140 to the proper operating strength and is then removed. 
[0044] The NMR signals picked up by the RF coil 152 are digitized by the transceiver module 150 and transferred 
to a memory module 1 60 which is also part of the system control 1 22. When the scan is completed and an entire array 

50 of data has been acquired in the memory module 1 60, an array processor 1 61 operates to Fourier transform the data 
into an array of image data. This image data is conveyed through the serial link 1 1 5 to the computer system 1 07 where 
it is stored in the disk memory 111 . In response to commands received from the operator console 100, this image data 
may be archived on the tape drive 112, or it may be further processed by the image processor 106 as will be described 
in more detail below and conveyed to the operator console 100 and presented on the video display 118. 

55 [0045] Referring particularly to Figs. 1 and 2, the transceiver 150 includes components which produce the RF exci- 
tation field B 1 through power amplifier 151 at a coil 152A and components which receive the resulting NMR signal 
induced in a coil 152B. As indicated above, the coils 152A and B may be separate as shown in Fig. 2, or they may be 
a single wholebody coil as shown in Fig. 1 . The base, or carrier, frequency of the RF excitation field is produced under 
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control of a frequency synthesizer 200 which receives a set of digital signals (CF) through the backplane 1 1 8 from the 
CPU module 119 and pulse generator module 121 . These digital signals indicate the frequency and phase of the RF 
carrier signal which is produced at an output 201 . The commanded RF carrier is applied to a modulator and up converter 
202 where its amplitude is modulated in response to a signal R(t) also received through the backplane 118 from the 

5 pulse generator module 121 . The signal R(t) defines the envelope, and therefore the bandwidth, of the RF excitation 
pulse to be produced. It is produced in the module 1 21 by sequentially reading out a series of stored digital values that 
represent the desired envelope. These stored digital values may, in turn, be changed from the operator console 100 
to enable any desired RF pulse envelope to be produced. The modulator and up converter 202 produces an RF pulse 
at the desired Larmor frequency at an output 205. 

w [0046] The magnitude of the RF excitation pulse output through line 205 is attenuated by an exciter attenuator circuit 
206 which receives a digital command, TA, from the backplane 118. The attenuated RF excitation pulses are applied 
to the power amplifier 151 that drives the RF coil 1 52A. For a more detailed description of this portion of the transceiver 
122, reference is made to U.S. Patent No. 4,952,877. 

[0047] Referring still to Fig. 1 and 2 the NMR signal produced by the subject is picked up by the receiver coil 152B 
15 and applied through the preamplifier 1 53 to the input of a receiver attenuator 207. The receiver attenuator 207 further 
amplifies the NMR signal and this is attenuated by an amount determined by a digital attenuation signal (RA) received 
from the backplane 118. The receive attenuator 207 is also turned on and off by a signal from the pulse generator 
module 121 such that it is not overloaded during RF excitation. 

[0048] The received NMR signal is at or around the Larmor frequency, which in the preferred embodiment is around 
20 63.86 MHz for 1.5 Tesla and 21.28 MHz for 0.5 Tesla. This high frequency signal is down converted in a two step 
process by a down converter 208 which first mixes the NMR signal with the carrier signal on line 201 and then mixes 
the resulting difference signal with the 2.5 MHz reference signal on line 204. The resulting down converted NMR signal 
on line 21 2 has a maximum bandwidth of 1 25 kHz and it is centered at a frequency of 1 87.5 kHz. The down converted 
NMR signal is applied to the input of an analog-to-digital (A/D) converter 209 which samples and digitizes the analog 
25 signal at a rate of 250 kHz. The output of the A/D converter 209 is applied to a digital detector and signal processor 
210 which produce 16-bit in-phase (I) values and 16-bit quadrature (Q) values corresponding to the received digital 
signal. The resulting stream of digitized I and Q values of the received NMR signal is output through backplane 1 1 8 to 
the memory module 160 where they are employed to reconstruct an image. 

[0049] To preserve the phase information contained in the received NMR signal, both the modulator and up converter 
30 202 in the exciter section and the down converter 208 in the receiver section are operated with common signals. More 
particularly, the carrier signal at the output 201 of the frequency synthesizer 200 and the 2.5 MHz reference signal at 
the output 204 of the reference frequency generator 203 are employed in both frequency conversion processes. Phase 
consistency is thus maintained and phase changes in the detected NMR signal accurately indicate phase changes 
produced by the excited spins. The 2.5 MHz reference signal as well as 5, 10 and 60 MHz reference signals are 
35 produced by the reference frequency generator 203 from a common 20 MHz master clock signal. The latter three 
reference signals are employed by the frequency synthesizer 200 to produce the carrier signal on output 201. For a 
more detailed description of the receiver, reference is made to U.S. Patent No. 4,992,736. 

[0050] Referring particularly to Fig. 3, a preferred embodiment of a pulse sequence which may be used to acquire 
NMR data according to the present invention is shown. Actually, two pulse sequences are shown, with the first being 
40 used to acquire NMR data for a synchronous spin motion image (A), and the second being used to acquire NMR data 
for a reference image (R). In the preferred embodiment, these two pulse sequences are alternated throughout the scan 
such that the corresponding views (i.e. phase encodings) in the data sets (A) and (R) are acquired at substantially the 
same moment in time. 

[0051] The pulse sequences are fundamentally a 2DFT pulse sequence using a gradient recalled echo. Transverse 
45 magnetization is produced by a selective 90° rf excitation pulse 300 which is produced in the presence of a slice select 
gradient (G z ) pulse 301 and followed by a rephasing gradient pulse 302. A phase encoding gradient (G y ) pulse 304 is 
then applied at an amplitude and polarity determined by the view number of the acquisition. A read gradient (G x ) is 
applied as a negative dephasing lobe 306, followed by a positive readout gradient pulse 307. An NMR echo signal 309 
is acquired 40 msecs. after the rf excitation pulse 300 during the readout pulse 307 to frequency encode the 256 
50 digitized samples. The pulse sequence is concluded with spoiler gradient pulses 312 and 313 along read and slice 
select axes, and a rephasing gradient pulse 311 is applied along the phase encoding axis (G y ). As is well known in the 
art, this rephasing pulse 31 1 has the same size and shape, but opposite polarity of the phase encoding pulse 304. The 
pair of pulse sequences are repeated 128 times with the phase encoding pulse 304 stepped through its successive 
values to acquire a 1 28 by 256 array of complex NMR signal samples that comprise the data set (A) and a 1 28 by 256 
55 array of complex NMR signal samples that comprise the reference data set (R). 

[0052] To practice the present invention an alternating magnetic field gradient is applied after the transverse mag- 
netization is produced and before the NMR signal is acquired. In the preferred embodiment illustrated in Fig. 3, the 
read gradient (G x ) is used for this function and is alternated in polarity to produce five bipolar, gradient waveforms 315. 
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The alternating gradient 315 has a frequency of 200 Hz and a duration of 25 msecs. At the same time, the pulse 
generator module 121 produces sync pulses as shown at 317, which are also at a frequency of 200 Hz and have a 
specific phase relationship with the alternating gradient pulses 31 5. As explained above, these sync pulses 31 7 activate 
the transducer 130 to apply an oscillating stress 319 to the patient which has the same frequency and phase relation- 
ship. To insure that the resulting waves have time to propagate throughout the field of view, the sync pulses 317 may 
be turned on well before the pulse sequence begins, as shown in Fig. 3. 

[0053] The phase of the NMR signal 309 A acquired during the first pulse sequence (A) is indicative of the movement 
of the spins. If the spins are stationary, the phase of the NMR signal is not altered by the alternating gradient pulses 
31 5, whereas spins moving along the read gradient axis (x) will accumulate a phase proportional to their velocity. Spins 
which move in synchronism and in phase with the alternating magnetic field gradient 215 will accumulate maximum 
phase of one polarity, and those which move in synchronism, but 1 80° out of phase with the alternating magnetic field 
gradient 21 5 will accumulate maximum phase of the opposite polarity. The phase of the acquired NMR signal 309 A is 
thus affected by system phase errors and random movement of spins along the x-axis, as well as the "synchronous" 
movement of spins along the x-axis. 

[0054] The reference pulse sequence is designed to measure the signal phase produced by sources other than 
synchronized spin movement. This is accomplished by repeating the identical pulse sequence, but without applying 
the oscillating stress 319. As a result, the phase of the acquired NMR signal 309 R will be affected by "static" system 
phase errors caused by field inhomogeneities and the like as well as the phase due to random spin movement along 
the x-axis. However, there will not be a phase component due to synchronous spin movement and the reference phase 
<|> R can, therefore, be subtracted from the phase <|> A to yield the phase (<|>) due solely to synchronous spin motion. 
[0055] The pulse sequence in Fig. 3 can be modified to measure synchronous spin movement along the other gradient 
axes (y and z). For example, the alternating magnetic field gradient pulses may be applied along the phase encoding 
axis (y) as indicated by dashed lines 321 , or they may be applied along the slice select axis (z) as indicated by dashed 
lines 322. Indeed, they may be applied simultaneously to two or three of the gradient field directions to "read" synchro- 
nous spin movements along any desired direction. 

[0056] The present invention may be implemented using most types of MR imaging pulse sequences. Gradient echo 
sequences can be readily modified to incorporate the alternating gradient as illustrated in the preferred embodiment. 
In some cases, however, the characteristics of a gradient echo sequence may not be ideal for a particular application 
of the technique. For example, some tissues (such as those with many interfaces between materials with dissimilar 
magnetic susceptibilities) may have a relatively short T2* relaxation time and therefore may not provide enough signal 
to obtain a noise-free image at the required echo delay time. In this setting, a spin echo implementation of the invention 
may be ideal, because for a given echo delay time TE, this pulse sequence is much less sensitive to susceptibility 
effects than a gradient echo sequence. When a spin echo pulse sequence is used, the alternating magnetic field 
gradient can be applied either before and/or after the 180° rf inversion pulse. However, if the alternating gradient is 
applied both before and after the rf inversion pulse, the phase of the alternating magnetic field gradient must be inverted 
180° after the rf inversion pulse in order to properly accumulate phase. 

[0057] In other applications, reduced acquisition time may be desirable. Fast spin echo and RARE sequences are 
rapid MRI sequences that acquire multiple views per TR cycle by applying different phase encoding gradients to each 
echo in a spin echo train. If 16 echoes, for instance, are acquired in each repetition of the sequence, then the total 
acquisition time for a complete image will be reduced by a factor of 16. One approach for modifying a fast spin echo 
sequence to implement the present invention is to insert the alternating gradient pulses between the initial 90° RF 
pulse and the first 1 80° RF refocussing pulse, followed by a similar but inverted set of gradient pulses. The first echo 
in the train might be at a TE of 40-60 msec, but the spacing between subsequent echoes could be as short as 12-15 
msec. 

[0058] Echo-planar imaging ("EPI") is another approach for high speed MR acquisition. In one version of this tech- 
nique, the spin echo created by standard 90° and 180° RF pulses is broken up into a series of 64-128 short gradient 
echoes by rapidly reversing the readout gradient. A different phase encoding is applied to each of the gradient echoes 
and therefore the acquired data from one shot of the sequence can, in principle, be used to reconstruct a complete 
image. One approach for modifying such an echo-planar sequence to implement the present invention is to insert 
alternating gradient pulses between the initial 90° RF pulse and the 180° RF refocussing pulse, followed by a similar 
but inverted set of alternating gradient pulses. Such a sequence permits phase images to be obtained in only a few 
seconds or less. 

[0059] The number of cycles of the alternating magnetic field gradient used in each pulse sequence depends on the 
strength of the applied gradient field, the frequency of the synchronous movement to be measured, and the TE time 
of the pulse sequence. The phase sensitivity of the pulse sequence to synchronous spin movement is proportional to 
the integrated product of alternating gradient field amplitude and the displacement over time. The sensitivity may be 
increased by increasing the amplitude of the gradient field pulses and by increasing the area under each pulse by 
making them as "square" as possible. The duration of each gradient pulse is limited by the desired synchronous fre- 
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quency, and hence more cycles of the alternating gradient waveform are required at higher frequencies to produce the 
same sensitivity as a lower frequency alternating gradient of the same amplitude and wave shape. 
[0060] Phase sensitivity to synchronous motion can also be increased by applying both the alternating gradient 
pulses 315 and the sync pulses 317 during the reference pulse sequence. However, when this is done the phase of 

5 the alternating magnetic field gradient 315 must be inverted 180° relative to the sync pulses 317 so that the sign of 
the accumulated phase is reversed. In addition, the magnetic field gradients should be flow compensated as described, 
for example, in U.S. patent No. 4,728,890 by Pattany et al. Such flow compensation removes any phase component 
due to random spin motion without significantly affecting the sensitivity to synchronous spin motion. Consequently, 
when the phase difference image is produced the phase accumulations due to synchronous spin motion add together, 

10 while phase accumulations due to other sources subtract and are thereby nulled. 

[0061 ] If the synchronous frequency to be measured exceeds the frequency at which the magnetic field gradient can 
be switched, a lower harmonic frequency may be used for the alternating gradient. As long as the time period of one 
cycle of the alternating gradient corresponds with an odd number of cycles of the synchronous spin motion, the phase 
of synchronous spin motion will accumulate. For example, the alternating gradient frequency may be 1/3, 1/5, 1/7, 1/9, 

15 etc. of the synchronous motion frequency. 

[0062] The oscillating stress may be applied by the transducer 1 30 in a number of ways. By starting the sync pulses 
31 7 well before the alternating magnetic field gradient 315 as shown in Fig. 3, the synchronous spin motion propagates 
throughout the field of view of the reconstructed image. This will image the steady-state conditions in the medium when 
the oscillating stress is applied. If the sync pulses 31 7 are turned off just before the alternating gradient 315 is applied, 

20 spins adjacent to the transducer 130 are moving with less amplitude or not at all during the phase accumulation time 
period. This may be desired, for example, when regions deep beneath the surface are of primary interest and large 
strain effects in the image near the transducer 1 30 can be suppressed. This is shown, for example, in Fig. 1 0a where 
the synchronous stress variations are depicted by the waveforms 336 and 338. The stress oscillations may be termi- 
nated one cycle short, as shown by waveform 336 when phase accumulation near the transducer 1 30 is to be reduced 

25 only slightly, or it may be terminated sooner as shown by waveform 338 if more phase accumulation is to be prevented. 
[0063] Referring particularly to Fig. 1 0b, the synchronous stress variations may only be applied for a moment in order 
to observe the mechanical impulse response or transient response of the gyromagnetic material. As indicated by wave- 
forms 340-343, the burst of stress oscillations may be applied at different times (t 0 , t 1 , t 2 , t 3 , etc.) before the application 
of the alternating magnetic field gradient pulses 315 and a two-dimensional strain image reconstructed at each timing. 

30 A three-dimensional data set is thus acquired in which one of the dimensions is the time at which the stress oscillations 
are applied. The response of the spins at any image pixel location in the field of view can then be plotted as a function 
of time. The same procedure may be used to measure the response of the medium when the oscillating stress is turned 
off, as shown by the waveforms 344-346 in Fig. 10b. 

[0064] In the preferred embodiment the frequency of the sync pulses 31 7 and the alternating magnetic field gradient 
35 315 are constant during each scan. Other information can be acquired, however, if the frequency of either is varied. 
For example, if the frequency of the alternating gradient 315 is swept through a series of frequencies during each pulse 
sequence, or a series of acquisitions of each view are acquired at different alternating gradient frequencies, the me- 
chanical response of the gyromagnetic material to a band of stress frequencies can be measured. Similarly, the sync 
pulses 31 7 can be swept through a series of frequencies during each pulse sequence or during a series of acquisitions 
40 of the same view to produce a band of synchronous spin movement frequencies. The same result may also be achieved 
by changing the shape of the applied oscillatory stress from substantially sinusoidal as depicted in Fig. 3 to some other 
shape such as a square wave, which contains higher frequency harmonics. 

[0065] Such multi-frequency excitations and detections may be required in some applications where reflections cause 
standing waves. One such application, for example, is imaging the head where standing waves are formed by reflections 
45 off the skull. By imaging with multiple frequencies and adding or averaging the resulting signals, the effects of such 
standing waves can be reduced. 

[0066] A scan using the pulse sequence of Fig. 3 is carried out under the direction of a program executed by the 
NMR system of Fig. 1 . Since many different mechanical properties can be imaged according to the present invention, 
many different programs may be stored and called up for use by the operator. Some of these programs are illustrated 
50 in the drawings by flow charts and described in detail below. However, it can be appreciated that many variations from 
these specific examples are possible. 

[0067] Referring particularly to Fig. 4, a scan may be performed according to the present invention to acquire NMR 
data from which ten images are reconstructed showing the propagation of strain through the gyromagnetic medium. 
The program for this scan is entered at 400 and the pulse sequence of Fig. 3 is downloaded to the pulse generator 
55 module 1 21 . The sync pulses 21 7 in this pulse sequence are timed to be in phase with the alternating motion encoding 
gradient 215 as indicated at process block 402. The pulse sequence is then performed the necessary number of times 
to acquire two complete NMR data sets as indicated at process block-404. As described above, one of these data sets 
(A) is phase encoded with the transducer 1 30 energized, and the other data set (R) is a reference that is acquired with 
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the transducer 130 deenergized. As will be described in detail below with reference to Figs. 6a-6c, a single phase 
image (<|>) is then reconstructed from these two NMR data sets (A) and (R) as indicated at process block 406. The 
intensity of each pixel in this phase image (<|)) is determined by the phase imparted to spins moving in synchronism 
with the applied oscillatory stress. This synchronous movement is a measure of the strain in the gyromagnetic material 
5 and the phase image (<|>) is, therefore, an image which shows by the brightness of its pixels the strain at corresponding 
locations (i.e. voxels) in the image plane. 

[0068] This phase image is a snapshot of the strain at one moment in the applied stress cycle, and to see the strain 
at other moments, the phase relationship of the sync pulses 217 is changed, as indicated at process block 408, and 
the acquisition and image reconstruction process is repeated. The scan continues to acquire NMR data and reconstruct 
10 strain images at successive increments (for example, 36°) of sync pulse phase. When an entire 360° cycle has been 
captured in successive strain images, the system branches at decision block 410. The set of strain images may then 
be viewed consecutively to "animate" the propagation of strain through the gyromagnetic medium as indicated at proc- 
ess block 412. 

[0069] The single or "animated" strain images obtained by this technique provide a view of the propagation of strain 

15 waves through the tissues imaged. Even without further processing, it is possible to qualitatively characterize tissue 
and to detect abnormalities by noting various features of the depicted waves such as wave-spacing, spatial phase, 
diffraction, and attenuation. The animated images are particularly useful for detecting small lesions such as tumors 
that differ in compliance from surrounding tissues, because the localized phase shift is more easily identified than in 
a single static strain image. A simple approach for integrating the information provided by these strain images with 

20 conventional modulus (anatomy) images provided by the same pulse sequence, is to multiply the corresponding gray 
scale values of the two types of images, thereby obtaining a composite image which provides both types of information. 
[0070] Although the primary focus of most medical applications of the invented imaging technique is to characterize 
tissue and to detect focal abnormalities such as tumors, the technique is also useful for characterizing acoustic systems 
perse. For instance, the technique may be used to characterize an array of acoustic energy sources that are intended 

25 for tissue ablation or lithotripsy. By exciting each acoustic driver individually and obtaining a strain image according to 
the present invention, the precise strain pattern associated with each transducer can be shown in situ, thereby showing 
the effects of diffraction and phase shifts caused by intervening tissue. These strain images could then be used to 
determine the proper phasing pattern to be applied to each of the elements in the acoustic drive array, so as to yield 
maximum strain at the target within the patient. As a final check, the tailored excitation so-produced could be applied 

30 to all of the drivers in the array at a low power level and a strain image could be obtained to confirm that the acoustic 
strain is maximized at the target tissue. The strain image could also be used to exclude the possibility of collateral 
damage to other locations when therapeutic power levels are applied to the acoustic drivers. 

[0071 ] While it is apparent that there are many medical applications of the invented imaging technique, it is also clear 
that the strain images provided by the technique are useful for a variety of industrial design and testing applications. 

35 For instance, it is a common requirement to assess the vibrational characteristics and stability of complex engine parts. 
Because it has not been possible to readily measure the internal strain caused by vibrational loading on such parts, it 
is often necessary to simulate such parameters using finite element analysis. The invented imaging method could be 
used to directly observe these phenomena. In this application, a model of the part would be made from a porous 
material that has a scaled approximation of mechanical properties of the real part. The pores of the model part would 

40 then be impregnated with a material such as oil, or water that would provide an MRI signal. Vibrational loading could 
then be applied, and the internal pattern of strain could be imaged using the invented method. As in medical applications, 
the objectives of non-medical applications include (1 ) characterization of material properties, (2) detection of focal flaws 
or abnormalities, (3) depicting the pattern of propagating waves in a transmitting medium, (4) characterizing an acous- 
tical driver or array of drivers. 

45 [0072] While strain images provide valuable information concerning the mechanical properties of the gyromagnetic 
media, the present invention enables one to measure and image many other mechanical properties as well. A scan 
which acquires the data necessary for such measurements may be carried out under the direction of a program illus- 
trated by the flow chart in Fig. 5. 

[0073] Referring particularly to Fig. 5, the scan begins at block 420 with the stress set to be applied along one axis. 
50 The direction of the alternating gradient is then set along the same axis at process block 422 and the phase of the 
sync pulses relative to the motion encoding gradient is set to zero as indicated at process block 424. The pulse sequence 
of Fig. 3 is then executed to acquire the two image data sets (A) and (R) at process block 426, as described above. A 
phase difference image (§) is then produced as indicated at process block 428, and the system then loops back at 
decision block 430. 

55 [0074] During the second pass, the sync pulse phase is offset 90° from the phase of the alternating gradient as 
indicated at process block 432, and another phase image (§) is acquired and reconstructed. As a result, after two 
passes as detected at decision block 430, two phase images indicative of strain S L and S Q are produced as defined 
above in respective equations (7) and (9). As discussed above, considerable information about the mechanical prop- 
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erties of the gyromagnetic media can be calculated from these two strain images alone. 

[0075] However, the scan may proceed as indicated at decision block 432 to acquire further data from which more 
exhaustive mechanical properties can be measured. More specifically, the direction of the alternating gradient is 
changed to align along an orthogonal axis as indicated at process block 434, and two more strain images are acquired 
and reconstructed. This may be repeated again as determined at decision block 432 with the alternating gradient 
switched in direction to a third orthogonal axis at process block 434. Thus, two strain images S L and S Q are produced 
with the alternating gradient directed along each of the three coordinate axes x, y and z to completely define the strain 
produced by the stress applied along one axis. 

[0076] As indicated at decision block 436, the scan may continue by altering the direction of the applied stress. In 
most applications, this is accomplished by switching the sync pulses to a different stress transducer 1 30 that is oriented 
in an orthogonal direction to the first transducer 130, although it is also possible to pause the scan at this point and 
instruct the operator to change the location of the transducer 130. In any case, the direction of the applied oscillatory 
stress is switched as indicated at process block 438 and the system loops back to process block 422 to repeat the 
above-described measurements. This may be repeated again for a third orthogonal stress axis, as determined at 
decision block 436. If the scan is exhaustively executed, the following eighteen images are acquired and reconstructed: 
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[0077] It should be apparent that not all of these images need be acquired during a scan. Instead, the operator 
selects the elastic property or properties to be imaged and the program is configured to acquire and reconstruct the 
necessary strain images. From these strain images the desired elastic property can be calculated at each pixel, as 
indicated at process block 440 and described above under the heading "General Description of The Invention". The 
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calculated values are then employed to modulate the intensity of its pixel location in an image. 
[0078] In the scans described above a phase image is produced from two acquired NMR data sets (A) and (R). 
Referring particularly to Figs. 6A-C, there are three methods which may be used for performing this function. As shown 
in Fig. 6A at process block 450, in the first method a two-dimensional Fourier transformation ("2DFT") is performed on 
5 each acquired data set (A) and (R). The phase (|V A is then calculated for each pixel location in the (A) image, and the 
phase (|)' R is calculated at each pixel of the (R) image as indicated at process block 452. 

V A = Tan' Q' A /r A 

10 

<>' fl = Tan 1 Q' R /I' R 

The phase difference (<])) at each image pixel is then calculated at process block 454 by subtracting the reference phase 

15 (|)' R from the synchronous motion encoded phase §' A . 

[0079] The second method is referred to in the art as the complex difference technique. Referring particularly to Fig. 
6B, the first step as indicated at process block 456 is to calculate the complex difference of the I and Q values of the 
corresponding sampled NMR signals in the two data sets (A) and (R). After these values are linearized, as indicated 
at process block 457, a two-dimensional Fourier transformation is performed on the resulting complex difference, as 

20 indicated at process block 458, and the phase <|> at each pixel in the resulting image is calculated as indicated at process 
block 460. 

[0080] The third method shown in Fig. 6C is very similar to the method of Fig. 6A. As indicated at process block 462, 
a two-dimensional Fourier transformation is performed on each NMR data set (A) and (R) and the complex difference 
between the resulting two images is calculated at process block 464. The resulting values are linearized as indicated 
25 at process block 465 and the phase of the resulting complex difference image is calculated at each pixel as indicated 
at process block 466. 

[0081 ] It should be apparent to those skilled in the art that many variations are possible in the particular method used 
to acquire the NMR data and reconstruct an image indicative of the spin phase accumulated as a result of synchronous 
spin motion. 

30 [0082] In the preferred embodiments described above, images are reconstructed in which pixel intensity is deter- 
mined directly by the phase of the synchronous motion encoded NMR signals. Other, less direct methods for detecting 
the synchronous motion are also possible. For example, if the frequency of the synchronous motion and alternating 
magnetic field gradient is increased, the spatial wavelength of the strain waves in normal tissues approaches the 
dimension of an imaged voxel. As a result, significant NMR signal phase dispersion occurs wit hin each voxel, and in 

35 a conventional image in which the brightness of each pixel is determined by the "modulus" (i.e. J? + Q 2 ) of the signal 
at each voxel, normal tissues will appear dark. On the other hand, the same strain waves passing through tumors and 
other less compliant tissues have a longer spatial wavelength and correspondingly less intravoxel phase dispersion. 
As a result, the signal intensity, or modulus, at these lesions will be much greater than the surrounding normal tissue 
and much brighter in the reconstructed image. 

40 [0083] Another less direct method for detecting the synchronous spin motion is to encode the resulting phase accu- 
mulation in longitudinal magnetization by tipping the phase encoded transverse magnetization back to the longitudinal 
axis. A pulse sequence which uses this synchronous phase encoded longitudinal magnetization will produce a modulus 
image that depicts the strain wave pattern in the subject. 

[0084] As indicated above, synchronous spin motion is produced using a transducer 130 that applies an oscillatory 
45 force to the gyromagnetic media. This force produces a corresponding stress in the gyromagnetic media which is 
propagated therethrough as strain in accordance with its elastic properties. This strain is the synchronous spin motion 
detected and imaged according to the present invention. 

[0085] Referring particularly to Figs. 7A and 7B, the preferred embodiment employs a transducer 130 suitable for 
applying an oscillatory force externally to a patient. It includes a coil of wire 510 wound on a bobbin 512, with the coil 

50 axis 514 directed perpendicular to the polarizing magnetic field B 0 . The coil 51 0 is 400 turns of 30 AWG. copper wire, 
and its leads 516 are connected directly to the amplifier 129 as described above. The bobbin 512 is mounted to a 
flexible arm 51 8 that is attached to, and cantilevered from a supporting block 520. When a current passes through the 
coil 510, the magnetic field which it produces interacts with the polarizing field B 0 . The bobbin 512 is thus twisted to 
bend the flexible arm 518 either upward or downward depending on the direction of current flow. By alternating the 

55 direction of current flow, the coil 51 0 twists back and forth to produce a corresponding alternating force which causes 
the flexible arm 51 8 to oscillate as indicated by the arrow 522. A strap 521 extends through an opening in the supporting 
block 520 and securely fastens it to the patient. 

[0086] The oscillatory motion of the arm 51 8 is coupled to the subject 530 by an applicator 524. In its simplest form 
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the applicator 524 is comprised of a tube 526 of the desired length which is secured at one of its ends to the arm 51 8. 
A pressure plate 528 is secured to the other end of the tube 526. The pressure plate 528 rests on the subject 530 to 
be imaged and its oscillatory movement produces a corresponding oscillatory compressive force that generates the 
desired synchronous spin motion. 

5 [0087] Other applicators 524 may be used with the transducer 130 to couple the oscillatory motion of the flexible 
arm 51 8 to the subject. For example, the applicator 524 may take the form of a probe which is inserted into an opening 
in the subject to rest against a particular structure (such as the prostate gland), or it may be inserted through tissue 
and its end anchored to an internal structure (such as the liver). Many shapes and sizes are possible. 
[0088] With the exception of the coil 510, all elements of the transducer 130 are constructed using non-magnetic 

10 materials. Thermoplastics such as acrylic may be used, for example, because such materials can be molded or ma- 
chined into many desired shapes and they are relatively benign to the tissues which they contact. 
[0089] For some applications it may be possible to produce the alternating gradient field pulses 315 with the same 
transducer 130 used to apply stress to the subject. In this case, the transducer 130 includes a coil that acts as a local 
gradient coil for surrounding tissues. The same alternating current applied to this coil to produce the alternating mag- 

15 netic field gradient pulses 31 5 also cause the coil to move in an oscillatory manner at the same frequency and phase. 
Such a coil might be used, for example, to detect tumors in the prostate gland. 

[0090] While application of a synchronous motion is preferably accomplished by the direct application of a force to 
the patient, indirect synchronous motion producing mechanisms are also possible. Electrical or other non-mechanical 
stimulation may be applied to the patient to produce synchronous motion in the tissues to be imaged. For example, 
20 synchronous brain function stimulation may produce synchronous motion at the cellular, or sub-cellular level. Synchro- 
nous electrical stimulation of selected muscles will, of course, produce synchronous motion in the muscles themselves, 
but also, detectable synchronous motion may be produced in the associated nerves and in tissues surrounding the 
muscles. 

[0091] Rather than applying a separate synchronous stress or synchronous stimulation, an endogenous periodic 
25 motion may also be used to practice the present invention. In this case the alternating motion encoding magnetic field 
gradient must be synchronized with the endogenous motion, and for best results, this requires the detection of a signal 
indicative of the frequency and phase of the endogenous motion. For example, the rapid movement of an abnormal 
heart valve may be detected with a phonocardiogram that senses the sound produced by the valve. The resulting 
electrical signal is used to synchronize the alternating magnetic field gradient. The resulting synchronous motion of 
30 the valve itself, as well as surrounding tissues and the blood flowing through the valve may be imaged. Similar results 
may be achieved with the heart using the ECG signal. 

[0092] Contrast agents are used in many imaging modalities to enhance image contrast in a diagnostically useful 
manner. Such a strategy is also possible with the present invention when the contrast agent operates to change the 
elastoviscous characteristics of the tissues of interest. For example, an osmotic agent which affects mechanical prop- 

35 erties in a manner that may be detected using the present invention may be used. A hyper-osmolar agent will draw 
water out of the target cells and make them more compliant, whereas a hypo-osmolar agent will add water to the target 
cells and make them rigid, or less compliant. This compliance difference can be imaged with the present invention and 
thus the affected cells can be contrasted against other cells that do not respond to the contrast agent. 
[0093] The present invention may also be used to enhance magnetic resonance angiograms ("MRA"). MRA relies 

40 on the movement of blood in the direction of one or more motion encoding magnetic field gradients to differentiate 
between stationary blood vessel walls and the blood which they carry. If the blood flow is slow, as in venous blood 
vessels, this contrast mechanism is less effective. The present invention may be used to produce an MRA by applying 
a compressive force at the synchronous frequency which sends compression waves through the vasculature. Since 
blood is an incompressible liquid, the synchronous compression wave will travel quite efficiently through the vessels 

45 as an oscillatory motion of the blood cells at the synchronous frequency. The alternating magnetic field gradient is 
"tuned" to this synchronous frequency and the resulting phase image will depict the vasculature. An example of this 
application of the invention will now be described. 

[0094] Referring particularly to Fig. 11, a cuff 550 is used in place of the transducer 130 in this embodiment of the 
invention. The cuff 550 receives the synchronous drive signal produced by the amplifier 129 (Fig. 1) at leads 551 and 

50 it constricts around the subject's ankle. The resulting synchronous compression waves in the venous structure and the 
arterial structure travel up the patient's leg 553 as indicated by arrows 555. The acquisition of a phase image as 
described above using the pulse sequence of Fig. 3 with a synchronous alternating magnetic field gradient 31 5 directed 
along the lengthwise dimension of the leg 553, may be used to produce an angiogram. To improve the SNR of the 
image additional phase images with the alternating magnetic field gradient directed along the other two orthogonal 

55 axes may also be acquired and summed together. 
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Claims 

1. A method of producing an image with an NMR imaging system which indicates a mechanical characteristic of 
tissue of a subject within the field of view of the NMR imaging system, the steps comprising: 

a) imparting (317, 129, 130) mechanical wave motion to spins located in the tissue 

by applying stress to the tissue with an alternating magnitude setting up a travelling mechanical wave in the 
tissue that causes cyclic strain in the stressed tissue; 

b) conducting (404, 426) a scan of the tissue with the NMR imaging system using a pulse sequence (Fig. 3) 
in which a magnetic field gradient (315) of alternating polarity is synchronized with the applied stress of alter- 
nating magnitude with a predetermined phase relationship and is applied to detect the cyclic wave motion of 
the spins, and 

c) reconstructing a phase image from NMR data (Fig. 6A -C) acquired during the scan (406, 428) by determining 
the phase of the NMR data at each pixel in a reconstructed image, thereby providing an image of the pattern 
of wave propagation in the tissue, and 

d) employing the phase information in the phase image to indicate the mechanical characteristic of the tissue 
at each pixel location of the phase image. 

2. The method as recited in claim 1 in which the mechanical wave motion is imparted to the spins by applying an 
oscillating stress (319) to the tissue. 

3. The method as recited in claim 1 or 2 to in which the mechanical wave motion is imparted to the spins by applying 
a stress (319) which varies periodically in magnitude. 

4. The method as recited in claim 1 , 2 or 3 in which the subject is a portion of a living animal. 

5. The method as recited in claim 2 in which the alternating magnetic field gradient (315) alternates a plurality of 
cycles before the acquisition of the NMR signal. 

6. The method as recited in claim 1 in which a second set of NMR data is acquired (404R, 426R) from the subject 
without imparting mechanical wave motion to the spins, and the phase indicated by the second set of NMR data 
is subtracted from the phase indicated by said reconstructed image (454, 456, 464). 

7. The method as recited in claim 1 in which a second phase image is produced from the subject by repeating (430) 
steps a), b) and c) with the relative phase between the mechanical wave motion and the alternating magnetic field 
gradient (315) being different than that used to produce said reconstructed image, and producing an image indic- 
ative of a mechanical characteristic of the subject (440) by combining said second phase image with said recon- 
structed image. 

8. The method as recited in claim 1 in which a second phase image is produced from the subject by repeating (432) 
steps a), b) and c) with the alternating magnetic field gradient (321, 322) being oriented in a different direction 
(434) than that used to produce said reconstructed image, and producing an image indicative of a mechanical 
characteristic of the subject (440) by combining the second phase image with said reconstructed image. 

9. The method as recited in claim 8 in which the respective alternating magnetic field gradients (31 5, 321 , 332) used 
to acquire NMR signals for the two images are orthogonal. 

10. The method as recited in claim 7 in which the mechanical characteristic (440) is one of a group of characteristics 
such as strain compliance, the velocity at which a strain wave is propagated in the subject, or the Young's modulus 
of the tissue. 

1 1 . The method as recited in claim 9 in which the mechanical characteristic (440) is one of a group of characteristics 
such as Poisson's ratio. 

12. The method as recited in claim 10 and 11, in which the bulk modulus or shear modulus are determined as said 
mechanical characteristic. 

13. An NMR system for producing an image which indicates a mechanical characteristic of tissue of a subject placed 
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within the field of view of the NMR system, which comprises: 

a transducer (130) arranged for coupling to the tissue for applying stress to spins located in the tissue; 
a magnet (140) for producing a polarizing magnetic field in the field of view; 
an RF transmitter (150) for producing transverse magnetization in the spins; 
an RF receiver (150) for receiving NMR signals produced by the spins; 

magnetic field gradient means (139) for producing a magnetic field in the field of view having a gradient in a 
first direction; 

a pulse generator (121) for producing a pulse sequence in which said transverse magnetization in the spins 
is produced, the magnetic field gradient means being adapted to produce a magnetic field gradient to detect 
the motion to of the spins, and the RF receiver (150) being adapted to receive NMR signals produced by the 
spins; and 

an image reconstructor (106) for reconstructing an image from the received NMR signals, characterized in 
that 

the transducer (1 30) is adapted to apply stress to the tissue with an alternating magnitude setting up a travelling 
mechanical wave in the tissue that causes cyclic strain in the stressed tissue; 

the magnetic field gradient means (139) is adapted to produce said magnetic field gradient (315) with an 
alternating polarity and being synchronized with the applied stress of alternating magnitude with a predeter- 
mined phase relationship for detecting the cyclic wave motion of the spins, and 

the image reconstructor (106) is adapted to determine the phase of the NMR data at each pixel in a recon- 
structed image, thereby providing an image of the pattern of wave propagation in the tissue, and is adapted 
to employ the phase information in the phase image to indicate the mechanical characteristic of the tissue at 
each pixel location of the phase image. 

14. The NMR system as recited in claim 13 in which the transducer (130) is adapted to impart the mechanical wave 
motion by applying a stress to the tissue which varies periodically in magnitude. 

15. The NMR system as recited in claims 13 or 14 in which the transducer (130) is adapted to impart the mechanical 
wave motion to the spins by applying an oscillating stress to a portion of the tissue. 

16. The NMR system as recited in any of the claims 13 to 15 in which the magnetic field gradient means (139) is 
adapted to alternate the polarity of the magnetic field gradient a plurality of cycles before receipt of the NMR signals. 

17. The NMR system as recited in claim 13 in which the pulse generator (121) is adapted to produce a second pulse 
sequence while the transducer (130) is turned off to receive second NMR signals, and the image reconstructor 
(1 06) is adapted to employ the phase of the received second NMR signals to indicate the mechanical characteristic 
of the subject. 

18. The NMR system as recited in claim 13 in which the pulse generator (121) is adapted to produce a second pulse 
sequence and to receive a second NMR signal in which the relative phase between the mechanical wave motion 
imparted by the transducer (130) and the alternating magnetic field gradient is different than the relative phase 
used for receiving said NMR signals, and the image reconstructor (106) is adapted to employ the phase of the 
received second NMR signals to indicate the mechanical characteristic of the subject. 

19. The NMR system as recited in claim 13 in which the pulse generator (121) is adapted to produce a second pulse 
sequence in which the magnetic field gradient means (139) is adapted to produce an alternating magnetic field 
gradient oriented in a different direction than that produced during said pulse sequence, and second NMR signals 
are received, and in which the image reconstructor is adapted to employ the phase of the received second NMR 
signals to indicate the mechanical characteristic of the subject. 

20. The NMR system as recited in claim 19 in which magnetic field gradient means is adapted to produce said alter- 
nating magnetic field gradient during the second pulse sequence with an orientation being orthogonal to the ori- 
entation of the alternating magnetic field gradient produced during said pulse sequence. 



Patentanspruche 

1. Verfahren zum Erzeugen eines Bildes miteinem NMR Abbildungssystem, welches eine mechanische Eigenschaft 
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vom Gewebe eines Objektes innerhalb des Betrachtungsfeldes des NMR Abbildungssystems anzeigt, aufweisend 
die Schritte: 

a) Beaufschlagen (317, 129, 130) von Spins, die im Gewebe liegen, mit mechanischer Wellenbewegung, 
5 indem das Gewebe mit einer alternierenden Amplitude beansprucht wird, welche im Gewebe eine sich aus- 

breitende mechanische Welle aufbaut, die im beanspruchten Gewebe eine zyklische Dehnung hervorruft; 

b) Durchfuhren (404, 426) einer Abtastung des Gewebes mit dem NMR Abbildungssystem unter Verwendung 
einer Impulssequenz (Fig. 3), in der ein Magnetfeldgradient (315) alternierender Polaritat mit der angewandten 
Beanspruchung alternierender Amplitude unter einer vorbestimmten Phasenbeziehung synch ronisiert ist und 

10 angelegt wird, urn die zyklische Wellenbewegung der Spins zu detektieren, und 

c) Rekonstruieren eines Phasenbildes aus NMR Daten (Fig. 6A - C), die wahrend der Abtastung (406, 428) 
erfaBt wurden, durch Ermitteln der Phase der NMR Daten bei jedem Pixel in einem rekonstruierten Bild, wo- 
durch ein Bild des Wellenausbreitungsmusters im Gewebe bereit gestellt wird, und 

d) Verwenden der-Phaseninformation im Phasenbild zur Anzeige der mechanischen Eigenschaften des Ge- 
15 webes bei jeder Pixelstelle im Phasenbild. 

2. Verfahren nach Anspruch 1, in welchem die mechanische Wellenbewegung den Spins dadurch auferlegt wird, 
daB das Gewebe einer oszillierenden mechanischen Beanspruchung (319) ausgesetzt wird. 

20 3. Verfahren nach Anspruch 1 oder 2, in welchem die mechanische Wellenbewegung den Spins dadurch auferlegt 
wird, daB eine Beanspruchung (319) angewandt wird, die periodisch in ihrer Amplitude variiert. 

4. Verfahren nach Anspruch 1 , 2 oder 3, in welchem der Gegenstand ein Teil eines lebendigen Wesens ist. 

25 5. Verfahren nach Anspruch 2, in welchem der alternierende Magnetfeldgradient (315) liber mehrere Zyklen vor der 
Erfassung des NMR Signals alterniert. 

6. Verfahren nach Anspruch 1 , in welchem ein zweiter Satz NMR Daten vom Objekt erfaBt wird (404R, 426R), ohne 
daB den Spins eine mechanische Wellenbewegung auferlegt wird, und die durch den zweiten Satz NMR Daten 

30 angezeigte Phase von der Phase subtrahiert wird, die durch das rekonstruierte Bild (454, 456, 464) angezeigtwird. 

7. Verfahren nach Anspruch 1 , in welchem ein zweites Phasenbild vom Objekt erzeugt wird, indem die Schritte a), 
b) und c) wiederholt werden (430) und dabei die relative Phase zwischen der mechanischen Wellenbewegung 
und dem alternierenden Magnetfeldgradienten (315) von derjenigen verschieden ist, die dazu verwendet wurde, 

35 das rekonstruierte Bild zu erzeugen, und in welchem ein Bild erzeugt wird, das eine mechanische Eigenschaft des 

Objekts (440) anzeigt, indem das zweite Phasenbild mit dem rekonstruierten Bild kombiniert wird. 

8. Verfahren nach Anspruch 1, in welchem ein zweites Phasenbild vom Objekt erzeugt wird, indem die Schritte a), 
b) und c) wiederholt werden (432), wobei der alternierende Magnetfeldgradient (321 , 322) in einer anderen Rich- 

40 tung (434) orientiert wird, als diejenige, die zur Erzeugung des rekonstruierten Bildes verwendet wurde, und in 

welchem ein Bild erzeugt wird, das eine mechanische Eigenschaft des Objektes (440) anzeigt, indem das zweite 
Phasenbild mit dem rekonstruierten Bild kombiniert wird. 

9. Verfahren nach Anspruch 8, in welchem die jeweiligen alternierenden Magnetfeldgradienten (315, 321, 332), die 
45 zur Erfassung von NMR Signalen fur die beiden Bilder verwendet werden, orthogonal sind. 

10. Verfahren nach Anspruch 7, in welchem die mechanische Eigenschaft (440) eine einer Gruppe von Eigenschaften 
ist wie Dehnung, Nachgiebigkeit, Geschwindigkeit, mit der eine Dehnungswelle sich im Objekt ausbreitet, oder 
der Elastizitatsmodul des Gewebes. 

50 

1 1 . Verfahren nach Anspruch 9, in welchem die mechanische Eigenschaft (440) eine einer Gruppe von Eigenschaften 
wie Poissonscher Beiwert ist. 

12. Verfahren nach Anspruch 10 und 11, in welchem der Kompressionsmodul oder Schermodul als die mechanische 
55 Eigenschaft ermittelt werden. 

13. NMR System zum Erzeugen eines Bildes, das eine mechanische Eigenschaft vom Gewebe eines Objektes anzeigt, 
das im Betrachtungsfeld des NMR Systems plaziert ist, aufweisend: 
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einen Wandler (1 30), der dazu ausgelegt ist, an das Gewebe angekoppelt zu werden, um Spins, die im Gewebe 
liegen, eine Belastung aufzuerlegen; 

einen Magneten (140) zum Erzeugen eines polarisierenden Magnetfeldes im Betrachtungsfeld; 

einen HF Transmitter (150) zum Erzeugen einer transversalen Magnetisierung in den Spins; 

einen HF Empfanger (150) zum Empfangen von NMR Signalen, die von den Spins hervorgerufen werden; 

eine Magnetfeldgradienteneinrichtung (1 39) zum Erzeugen eines Magnetfeldes im Betrachtungsfeld mit einem 

Gradienten in einer ersten Richtung; 

einen Impulsgenerator (121) zum Erzeugen einer Impulssequenz, in der die transversale Magnetisierung in 
den Spins hervorgerufen wird, wobei die Magnetfeldgradienteneinrichtung dazu ausgelegt ist, einen Magnet- 
feldgradienten zur Detektion der Bewegung der Spins zu erzeugen, und der HF Empfanger (150) dazu aus- 
gelegt ist, NMR Signale, die von den Spins hervorgerufen werden, zu empfangen; und 
eine Bildrekonstruktionseinrichtung (106) zum Rekonstruieren eines Bildes aus den empfangenen NMR Si- 
gnalen, dadurch gekennzeichnet, daft 

der Wandler (130) dazu ausgelegt ist, das Gewebe mit einer Beanspruchung alternierender Amplitude zu 
beaufschlagen, die im Gewebe eine sich ausbreitende mechanische Welle aufbaut, welche im beanspruchten 
Gewebe eine zyklische Dehnung hervorruft; 

die Magnetfeldgradienteneinrichtung (139) dazu ausgelegt ist, den Magnetfeldgradienten (315) mit einer al- 
ternierenden Polaritat zu erzeugen und so, daB er mit der beaufschlagten Beanspruchung alternierender Am- 
plitude unter einer vorbestimmten Phasenbeziehung zur Detektion der zyklischen Wellenbewegung der Spins 
synchronisiert ist, und 

die Bildrekonstruktionseinrichtung (106) dazu ausgelegt ist, die Phase der NMR Daten bei jedem Pixel in 
einem rekonstruierten Bild zu bestimmen, wodurch ein Bild des Wellenausbreitungsmusters im Gewebe bereit 
gestellt wird, und dazu ausgelegt ist, die Phaseninformation im Phasenbild dazu auszunutzen, die mechani- 
sche Eigenschaft des Gewebes bei jeder Pixelstelle im Phasenbild anzuzeigen. 

14. NMR System nach Anspruch 13, in welchem der Wandler (130) dazu ausgelegt ist, die mechanische Wellenbe- 
wegung dadurch aufzuerlegen, daB er das Gewebe einer Beanspruchung unterzieht, die periodisch in der Ampli- 
tude variiert. 

15. NMR System nach Anspruch 13 oder 14, in welchem der Wandler (130) dazu ausgelegt ist, die mechanische 
Wellenbewegung den Spins dadurch aufzuerlegen, daB er auf einen Teil des Gewebes eine oszillierende Bean- 
spruchung ausubt. 

16. NMR System nach einem der Anspruche 13 bis 15, in welchem die Magnetfeldgradienteneinrichtung (139) dazu 
ausgelegt ist, die Polaritat des Magnetfeldgradienten mehrere Zyklen vor Empfang der NMR Signale zu alternieren. 

17. NMR System nach Anspruch 13, in welchem der Impulsgenerator (121) dazu ausgelegt ist, eine zweite Impuls- 
sequenz zu erzeugen, wahrend der Wandler (130) abgeschaltet ist, um zweite NMR Signale zu empfangen, und 
die Bildrekonstruktionseinrichtung (106) dazu ausgelegt ist, die Phase der empfangenen zweiten NMR Signale 
dazu heranzuziehen, die mechanische Eigenschaft des Objektes anzuzeigen. 

18. NMR System nach Anspruch 13, in welchem der Impulsgenerator (121) dazu ausgelegt ist, eine zweite Impuls- 
sequenz zu erzeugen und ein zweites NMR Signal zu empfangen, in welchem die relative Phase zwischen der 
mechanischen Wellenbewegung, die vom Wandler (130) auferlegt wird, und dem alternierenden Magnetfeldgra- 
dienten sich von der relativen Phase unterscheidet, die zum Empfang dieser NMR Signale verwendet wurde, und 
in welchem die Bildrekonstruktionseinrichtung (1 06) dazu ausgelegt ist, die Phase der empfangenen NMR Signale 
zur Anzeige der mechanischen Eigenschaft des Objektes zu verwenden. 

19. NMR System nach Anspruch 13, in welchem der Impulsgenerator (121) dazu ausgelegt ist, eine zweite Impuls- 
sequenz zu erzeugen, in der die Magnetfeldgradienteneinrichtung (1 39) dazu eingerichtet ist, einen alternierenden 
Magnetfeldgradienten zu erzeugen, der in einer anderen Richtung als der Richtung orientiert ist, die wahrend 
dieser Impulssequenz erzeugt wurde, und in welchem zweite NMR Signale empfangen werden, und in welchem 
die Bildrekonstruktionseinrichtung dazu ausgelegt ist, die Phase der zweiten empfangenen NMR Signale zur An- 
zeige der mechanischen Eigenschaft des Objektes zu verwenden. 

20. NMR System nach Anspruch 19, in welchem die Magnetfeldgradienteneinrichtung dazu ausgelegt ist, den alter- 
nierenden Magnetfeldgradienten wahrend der zweiten Impulssequenz mit einer Orientierung zu erzeugen, die zur 
Orientierung des alternierenden Magnetfeldgradienten orthogonal ist, der wahrend dieser Impulssequenz erzeugt 
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wurde. 



Revendications 

1. Procede de production d'une image avec un systeme d'imagerie par resonance magnetique nucleaire (RMN) qui 
indique une caracteristique mecanique du tissu d'un sujet dans le champ d'observation du systeme d'imagerie 
RMN, dont les etapes comprennent : 

a) la communication (31 7, 1 29, 1 30) d'un mouvement mecanique ondulatoire a des spins situes dans le tissu 
en appliquant une contrainte au tissu avec une amplitude alternative, etablissant une onde mecanique pro- 
gressive dans le tissu qui provoque une deformation cyclique dans le tissu sous contrainte ; 

b) la realisation (404, 426) d'un balayage du tissu avec le systeme d'imagerie RMN en utilisant une sequence 
d'impulsions (figure 3) dans laquelle un gradient de champ magnetique (315) de polarite alternative est syn- 
chronise avec la contrainte appliquee d'amplitude alternative avec un rapport de phase predetermine et est 
applique afin de detecter le mouvement cyclique ondulatoire des spins ; 

c) la reconstruction d'une image de phase a partir des donnees RMN (figures 6A-C) acquises pendant le 
balayage (406, 428) en determinant la phase des donnees RMN a chaque pixel dans une image reconstruite, 
afin de fournir ainsi une image du schema de propagation d'onde dans le tissu ; et 

d) I'emploi des informations de phase dans I'image de phase afin d'indiquer la caracteristique mecanique du 
tissu a chaque position de pixel de I'image de phase. 

2. Procede selon la revendication 1, dans lequel le mouvement mecanique ondulatoire est communique aux spins 
en appliquant une contrainte oscillatoire (319) au tissu. 

3. Procede selon la revendication 1 ou 2, dans lequel le mouvement mecanique ondulatoire est communique aux 
spins en appliquant une contrainte (319) qui varie periodiquement en amplitude. 

4. Procede selon la revendication 1 , 2 ou 3, dans lequel le sujet est une partie d'un animal vivant. 

5. Procede selon la revendication 2, dans lequel le gradient de champ magnetique alternatif (315) effectue des al- 
ternances pendant plusieurs cycles avant I'acquisition du signal RMN. 

6. Procede selon la revendication 1 , dans lequel un deuxieme jeu de donnees RMN est acquis (404R, 426R) depuis 
le sujet sans communiquer de mouvement mecanique ondulatoire aux spins, et la phase indiquee par le deuxieme 
jeu de donnees RMN est soustraite de la phase indiquee par ladite image reconstruite (454, 456, 464). 

7. Procede selon la revendication 1, dans lequel une deuxieme image de phase est produite depuis le sujet en 
repetant (430) les etapes a), b) et c) avec la phase relative entre le mouvement mecanique ondulatoire et le 
gradient de champ magnetique alternatif (315) qui est differente de celle utilisee pour produire ladite image re- 
construite, et en produisant une image indicative d'une caracteristique mecanique du sujet (440) en combinant la 
deuxieme image de phase avec ladite image reconstruite. 

8. Procede selon la revendication 1, dans lequel une deuxieme image de phase est produite depuis le sujet en 
repetant (432) les etapes a), b) et c) avec le gradient de champ magnetique alternatif (321 , 322) qui est oriente 
dans une direction differente (434) de celle utilisee pour produire ladite image reconstruite, et en produisant une 
image indicative d'une caracteristique mecanique du sujet (440) en combinant la deuxieme image de phase avec 
ladite image reconstruite. 

9. Procede selon la revendication 8, dans lequel les gradients de champ magnetique alternatif respectifs (315, 321 , 
332) utilises pour saisir des signaux RMN pour les deux images sont orthogonaux. 

10. Procede selon la revendication 7, dans lequel la caracteristique mecanique (440) appartient a un groupe de ca- 
racteristiques telles que la deformation, I'elasticite, la vitesse a laquelle une onde de deformation se propage dans 
le sujet, ou le module de Young du tissu. 

11. Procede selon la revendication 9, dans lequel la caracteristique mecanique (440) appartient a un groupe de ca- 
racteristiques telles que le coefficient de Poisson. 
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12. Procede selon la revendication 10 ou 11, dans lequel le module de compressibility ou le module de cisaillement 
sont determines en tant que ladite caracteristique mecanique. 

13. Systeme RMN pour produire une image qui indique une caracteristique mecanique du tissu d'un sujet place dans 
5 le champ d'observation du systeme RMN, qui comprend : 

un transducteur (1 30) dispose pour etre couple au tissu pour appliquer une contrainte a des spins situes dans 
le tissu ; 

un aimant (140) pour produire un champ magnetique polarisant dans le champ d'observation ; 
10 - un emetteur RF (150) pour produire une aimantation transversale dans les spins ; 

un recepteur RF (150) pour recevoir des signaux RMN produits par les spins ; 

des moyens a gradient de champ magnetique (139) pour produire un champ magnetique dans le champ 
d'observation ayant un gradient dans une premiere direction ; 

un generateur d'impulsions (121) pour produire une sequence d'impulsions dans laquelle ladite aimantation 
15 transversale dans les spins est produite, les moyens a gradient de champ magnetique etant adaptes pour 

produire un gradient de champ magnetique afin de detecter le mouvement des spins, et le recepteur RF (1 50) 
etant adapte pour recevoir des signaux RMN produits par les spins ; et 

un reconstructeur d'image (1 06) pour reconstruire une image a partir des signaux RMN recus, 

20 caracterise en ce que : 

le transducteur (130) est adapte pour appliquer une contrainte au tissu avec une amplitude alternative eta- 
blissant une onde mecanique progressive dans le tissu qui provoque une deformation cyclique dans le tissu 
sous contrainte ; 

25 - les moyens a gradient de champ magnetique (139) sont adaptes pour produire ledit gradient de champ ma- 

gnetique (315) avec une polarite alternative et qui est synchronise avec la contrainte appliquee d'amplitude 
alternative avec un rapport de phase predetermine pour detecter le mouvement cyclique ondulatoire des 
spins ; et 

le reconstructeur d'image (1 06) est adapte pour determiner la phase des donnees RMN a chaque pixel dans 
30 une image reconstruite, afin de fournir ainsi une image du schema de propagation d'onde dans le tissu, et est 

adapte pour employer les informations de phase dans I'image de phase afin d'indiquer la caracteristique me- 
canique du tissu a chaque position de pixel de I'image de phase. 

14. Systeme RMN selon la revendication 13, dans lequel le transducteur (130) est adapte pour communiquer le mou- 
35 vement mecanique ondulatoire en appliquant une contrainte au tissu qui varie periodiquement d'amplitude. 

15. Systeme RMN selon la revendication 13 ou 14, dans lequel le transducteur (130) est adapte pour communiquer 
le mouvement mecanique ondulatoire aux spins en appliquant une contrainte oscillatoire a une partie du tissu. 

40 16. Systeme RMN selon I'une quelconque des revendications 13 a 15, dans lequel les moyens a gradient de champ 
magnetique (1 39) sont adaptes pour faire alterner la polarite du gradient de champ magnetique pendant plusieurs 
cycles avant la reception des signaux RMN. 

17. Systeme RMN selon la revendication 13, dans lequel le generateur d'impulsions (121) est adapte pour produire 
45 une deuxieme sequence d'impulsions tandis que le transducteur (130) est eteint afin de recevoir des deuxiemes 

signaux RMN, et le reconstructeur d'image (1 06) est adapte pour employer la phase des deuxiemes signaux RMN 
recus afin d'indiquer la caracteristique mecanique du sujet. 

18. Systeme RMN selon la revendication 13, dans lequel le generateur d'impulsions (121) est adapte pour produire 
50 une deuxieme sequence d'impulsions et pour recevoir un deuxieme signal RMN dans lequel la phase relative 

entre le mouvement mecanique ondulatoire communique par le transducteur (130) et le gradient de champ ma- 
gnetique alternatif estdifferente de la phase relative utilisee pour recevoir lesdits signaux RMN, et le reconstructeur 
d'image (106) est adapte pour employer la phase des deuxiemes signaux RMN recus afin d'indiquer la caracte- 
ristique mecanique du sujet. 

55 

19. Systeme RMN selon la revendication 13, dans lequel le generateur d'impulsions (121) est adapte pour produire 
une deuxieme sequence d'impulsions dans laquelle les moyens agradientde champ magnetique (139) sont adap- 
tes pour produire un gradient de champ magnetique alternatif oriente dans une direction differente de celle produite 
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pendant ladite sequence d'impulsions et des deuxiemes signaux RMN sont recus, et dans lequel le reconstructeur 
d'image est adapte pour employer la phase des deuxiemes signaux RMN recus afin d'indiquer la caracteristique 
mecanique du sujet 

20. Systeme RMN selon la revendication 19, dans lequel les moyens a gradient de champ magnetique sont adaptes 
pour produire ledit gradient de champ magnetique alternatif pendant la deuxieme sequence d'impulsions avec une 
orientation qui est orthogonale a I'orientation du gradient de champ magnetique alternatif produit pendant ladite 
sequence d'impulsions. 
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